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Preface 



The last decade witnessed a radical change in cosmology - the science of the universe. 
Cosmology now becomes an observation dependent science, like other branches of physics. 
This is brought about by the high precision observation techniques developed over the last 
few years. The great surprise that results from the high precision data is the inference that the 
universe is undergoing an accelerated expansion at present defying all intuitions. The search 
for the matter responsible for this unexpected behaviour of the universe provides one of the 
greatest excitements in contemporary theoretical physics. 

The present thesis is a collection of five papers based on my research on the problem of 
'dark energy', the driver of alleged present acceleration of the universe. All the five papers are 
published in international journals. 

The thesis is divided into five chapters. The first chapter is an introduction where the 
problem is defined and a survey of the work already in the literature has been made. A brief 
outline of the present work is also included in the introduction. 

The next four chapters include the actual work done. The reprints of the published papers 
are presented as the chapters or sections thereof. Only the last chapter contains a slightly 
improved version of the published paper. 



Sudipta Das 



Chapter 1 
Introduction 



1 



Introduction 2 

The universe consists of everything that we can see through our best gadgets, have seen in 
the past and expect to see in any intelhgible future. Cosmology deals with the physics of this 
most exhaustive collection of objects as a whole. Cosmology is concerned with the formation, 
the evolution, the future of the universe made of billions of galaxies spread over billions of 
light years. 

Cosmology really came of age as a science after the advent of general relativity in 1915 
which made possible a systematic theoretical modelling of the universe and Bubble's observa- 
tion in 1929 that the universe is in fact evolving and hence interesting as a physical science. 
Like all other branches of physics, cosmology is also an observational science, but until very 
recently, the data available had been limited and that too had been plagued by the lack of 
precision. The dramatic change in the scenario started with the COSMIC BACKGROUND 
EXPLORER (COBE) |T1 [21 [Sj H] and over the past decade there had been an explosion of 
really high precision data, such as those from WILKINSON MICROWAVE ANISOTROPY 
PROBE (WMAP) [H O [71 El [9] . One inherent problem in cosmology, that one cannot repeat 
experiments, will remain for ever, one cannot ask the universe to evolve afresh with various 
initial conditions. But the data available on the existing system holds the key to dictate the 
direction of research in this branch. 

Cosmology as a subject is as vast as the universe, and application of every branch of 
physics is indeed warranted. The recent advances in observational cosmology has led to many 
exciting discoveries and possibilities regarding the evolution of the universe, but arguably 
the most exciting and puzzling amongst them is that the universe now is expanding with 
an acceleration defying the properties of the known matter content of the universe. For a 
systematic and lucid review of the observational results and their interpretations, we refer to 
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the work by L. Perivolaropoulos [10]. The present thesis endeavours to look at this problem 
from a very narrow angle. 

As already mentioned, this acceleration is counter intuitive as gravity, which is the deciding 
interaction in the governance of the dynamics of the universe as a whole, is always attractive. 
Gravity is by far the weakest amongst the four basic interactions of nature, but the strong 
and weak interactions are short range and have hardly anything to do with the dynamics of 
the universe at a large scale, while electromagnetic interaction has almost no impact in view 
of the charge neutrality of the universe. So the galaxies should attract each other, and even 
if the universe expands as suggested by the observations as well as the standard models of 
cosmology, the rate of expansion should be decreasing. 



1.1 Standard Cosmological Model : 

For an observationally realistic and logically viable model of the universe, one makes the fol- 
lowing assumptions: 

(1) General Relativity (GR) correctly describes gravity. 

(2) Cosmological principle is valid, i.e, universe on a large scale (> 10^ light years) is spatially 
homogeneous and isotropic. The typical size of a galaxy is ~ 10^ light years, an order of 
magnitude beyond which there is no preferred position or direction in the universe. 

(3) Hydrodynamic approximation : According to which the basic building blocks of the uni- 
verse are galaxies and their distribution can be considered as a fluid distribution. 

The most general metric satisfying the homogeneity and isotropy of the universe is given 
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by 

ds^ ^ dt^ - a\t)[- —+r^d9^ + r\m'^9d(f)^] , (1.1) 

1 — kr^ 

where t is the time and r, 9, are space co-ordinates, a{t) is the scale factor of the universe 
which gives the expansion history of the universe and k is called the curvature index. This 
metric is called the Friedmann-Robertson- Walker (FRW) metric. 

With the hydrodynamic approximation that the matter distribution of the universe can be 
approximated to be a perfect fluid, the energy momentum tensor for a perfect fluid distribution 
is taken as 

Tfj,u = (p + p)v^v^ - pg^y , (1.2) 

where t'^'s are the components of the fluid velocity vector, p is the energy density and p is the 
isotropic pressure of the perfect fluid. 
With this input, the Einstein equations 

= SttGT^, (1.3) 

lead to the differential field equations 

3^ + 3-1 = SttGa (1-4) 

a d"^ + k ^ , ^ 

2- + ^ = -SttGp, 1.5 

where a dot denotes differentiation w.r.t. the cosmic time t. 

Also we obtain a third equation, called the matter conservation equation, of the form 

p + 3-(p + ij) = 0, (1.6) 
a 

which is not an independent equation but can be derived from the two field equations or from 
the Bianchi identities. It deserves mention that the curvature index k can not be determined 
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from the field equations, and is rather put in by hand and can take values 0, +1 and -1 which 
correspond to flat, closed and open universes respectively. 

This system of equations can not be solved completely as there are three unknowns a, p and 
p and only two independent equations. So there is the need of a third equation for solving 
the system completely which is provided by the equation of state connecting the density and 
pressure of the cosmic fluid as 



Normally this equation of state is taken as that of a barotropic fluid, i.e, w is a constant. 
Thus, at present when the universe is expected to be matter dominated where there is no 
pressure (p = 0), w = 0. Similarly, at very early epoch, when the universe was very hot and 



be solved for and thus one can get an idea about the evolution of the universe. 

To relate this FRW model with the observations, some useful parameters are defined, 
(i) Hubble constant H is defined as 



which is an observable parameter and gives the expansion rate of the universe. This parameter 
has the dimension of (time)~^. 




(1.7) 



was dominated by radiation, then p = |p, i.e, w 



So, from these set of equations, a{t) can 




(1.8) 



(ii) Although the universe is expanding, it is dominated by the gravitational interaction 
which gives rise to an attractive force. So this expansion is expected to be decelerated. To 
express this deceleration, a dimensionless deceleration parameter q is defined as 
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g > indicates that ^ is negative, i.e, universe is decelerating and g < indicates that ^ 
is positive, i.e, universe is accelerating. 



(iii) A jerk parameter r is defined as 

d/a 

r = ^, 1.10 

which is a measure of the rate of change of q. As present observations facilitate the study of 
the evolution of the deceleration parameter q, this jerk parameter has become useful. Sahni et 
al [11] introduced a pair of ''statefinder parameters" to characterize the quintessence models, 
particularly the interacting models. In addition to r, the other parameter s is defined as 

r — 1 

^=Tw TT- (1-11) 

3(g- 2) 

(iv) The density of the universe is expressed in a dimensionless form by defining a density 
parameter given by 

fi = -^, (1.12) 

Pc 

where pc = |^ is called the critical density or closure density of the universe. 
In terms of these parameters, equations (11.41) and (II. 5p can be expressed as 



k 



H\l -2q) + — = -SttGp . (1.14) 

Now, if one considers a matter dominated universe where p = 0, equations (11.131) and (I1.14p 
leads to three possibilities : 

(a) When k = —1, q < ^ and fl < 1, i.e, p < pc- This means the matter density being less 
than the critical density, the universe will go on expanding for ever. Such models are called 
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i 

Figure 1.1: Evolution of the scale factor with time for open, fiat and closed models 
open universe models. 

(b) When k — +1, q > | and Q > 1, i.e, p > pc- So, the matter density being higher 
than the critical density, the universe will expand upto some maximum volume and then due 
to gravitational attraction it will re-collapse. Such models of the universe are called closed 
universe models. 

(c) When k — 0, q = ^ and D, — 1, i.e, p = Pc, one has the limiting case between the above 
two. This is called a flat model as the space section has zero curvature. 

In 1929, Hubble made a remarkable discovery regarding the motion of the galaxies. He 
observed that the galaxies are moving away from each other, and the velocity of separation v 
between two galaxies is proportional to the distance D between them, i.e, 

V oo D , 

or, v^HD , (1.15) 
where H — - is the Hubble constant defined earher. This means the galaxies were closer 
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together earlier. So, by tracking back one can arrive at a time when all the matter were con- 
centrated at a single point. At that instant universe had a zero volume and an infinite density. 
That epoch, when a{t) = and if — > oo corresponds to some violent activity and is given the 
name Big Bang Singularity. The existence of singularity in a theory is an unwanted feature 
as laws of physics break down and naturally the features cannot be explained. But still the 
Big Bang model enjoys the status of a preferred theory as it has its own success stories : 

(i) Can predict He abundance : One of the fundamental problems of cosmology is to explain 
the primary creation of matter and to understand the observed abundances of different ele- 
ments. At the beginning, when all the matter-energy of the universe was concentrated in a 
tiny volume, the spectrum of particles that we see today was surely absent. Although the 
singular stage is definitely out of our purview of explanation, the Big Bang theory can in fact 
trace the history of the universe when its size was around 10"'^^ cm. This is clearly much 
shorter than the de Broglie wavelength of most of the particles that we see today. Following 
Gamow's seminal work [12], one can explain the nucleosynthesis process that took place during 
the radiation dominated era. Particularly the abundance of lighter elements like Helium are 
well explained in Big Bang theory and the theoretical prediction has a close semblance with 
the observational results. 

(ii) Could predict the relic thermal radiation at ~ 2.7 °K. This Cosmic Microwave Background 
Radiation (CMBR) was detected later in 1965 by Penzias and Wilson \T3\. This detection 
confirmed the assumption of isotropy of the universe and is considered to be the greatest 
triumph of the Big Bang theory and is arguably the strongest pillar of modern cosmology. 

(iii) Age of the universe : The Friedmann models could provide a formula for the age of the 
universe which goes as 
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t - 



^ for a flat matter dominated universe. 



Thus knowing the value of Hq ( the subscript '0' indicates the present time ), to can be easily 
calculated. Following this method, the standard Big Bang theory predicts the correct order 
of magnitude of the age of the universe ~ 1.5 X 10^'' years. 

(iv) Formation of galaxies : Although the universe is homogeneous at a large scale, there are 
clumps ( galaxies and clusters of galaxies ) around us. So at a smaller scale, there should 
be inhomogeneity. The universe starts homogeneous, and still looks homogeneous at a scale 
more than 10^ light years, but at a smaller scale must have inhomogeneities, i.e, structures 
like galaxies. In the purview of standard Big Bang cosmology, this can also be explained as 
it has been shown that if there is some kind of perturbation, it can indeed give rise to some 
growing mode so that galaxies are formed. 

1.2 Problems of Standard Cosmological Model : 

Standard Big Bang cosmology has its share of problems too. A few of them are : 

(1) Horizon problem : Given the present size of the universe (~ 10^° light years), it is quite 
possible to consider two far separated points such that there is no causal connection between 
these two points, i.e, their light cones never intersect even if one traces back to the last 
scattering surface ( LSS ) when matter and radiation decoupled and the universe became 
transparent. But even then the two points carry the same information at present as the 
universe is homogeneous and isotropic. This is known as the ^horizon problem\ 

(2) Flatness problem : From equation fll.lSp . one can arrive at a relation 



n-1 



k 
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If it is considered that the initial conditions including the density parameter Q were set during 
the GUT epoch when temperature of the universe was ~ 10^^ GeV, then 

n = l±S, where S < 10"^° . 

This means the departure from the f2 = 1 value has to be very small. Any relaxation from 
this fine-tuning would have led to a much higher (or lower) value of Q at present which is not 
obtained observationally. 

This fine-tuned value of ^ 1 leads to a = model, i.e, a spatially flat model of the 
universe. Without this extreme fine tuning, the universe would have either collapsed back 
within a time scale of 10"^^ s (in a closed model) or would have expanded at a much higher 
rate (in an open model) than observed at present. Standard Big Bang model can not explain 
why fl is so closely tuned to 1 and the problem is termed the 'flatness problem' or equivalently 
the 'fine tuning problem\ 

(3) The monopole problem : Gauge field theories suggest that whenever there is any symmetry 
breaking, inevitably some particles are created which have the characteristics of magnetic 
monopoles. So, it is expected that during the phase transition of the universe, some monopoles 
must have been created which being highly stable particles, should have been observed at 
present epoch also. But, in practice monopoles are not observed. This is known as the 
'monopole problem\ 



1.3 The Inflationary Paradigm : 

The solution to this problem was suggested by Alan Guth in 1981 [H] by introducing the 
so called inflationary model of the universe. In this model, Guth suggested that during a 
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very early epoch the universe had a very rapid phase of accelerated expansion having quite 
a number of e-foldings of the volume in a short span of time. At that time universe was 
dominated by vacuum energy and the equation of state was of the form 

Pvac ~l~ Pvac ; 

where pyac is the vacuum energy density and p^ac is the corresponding pressure. 
Now, from the conservation equation (11.61) . one obtains p^ac = —Pvac =constant. Thus A = 
SiiGpvac serves as an effective cosmological constant. Then Einstein field equations (11. 4p and 
(ll.Sp can be written as 

3^ + 34 = A, (1-16) 

a + k , , , 

2- + — ^ = A. 1.17 

Equations (11.161) and (I1.17P can be combined to yield 

a = . 

For a positive A, a is thus positive and the universe has an accelerated expansion which is 
termed as the inflationary scenario. The deceleration parameter g = — 1| is negative definite. 
From equation (11.131) one can write 

VL = 2qH{VL~l) , (1.18) 

which clearly indicates that at 17 = 1, one has 17 = 0, i.e, 17 = 1 is the stable solution. For 
a negative value of g, 17 < for 17 > 1 and 17 > for 17 < 1 . So the density parameter 
decreases or increases to the stable value of 17 = 1 when it is greater or less respectively than 
17 = 1. Hence the contribution from the spatial curvature, 17^ = —-^^Jj^, is washed out in view 
of equation (ll.lSp . Hence, flatness problem was solved. 
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The horizon problem was also solved in the following way. Two points, which were causally 
connected during a very early epoch, might have fallen so much apart during the inflationary 
expansion that at this epoch their past light cones do not have any intersection even if they 
are extended back to the last scattering surface. 

Inflationary models could provide solution to the monopole problem also by considering 
that the monopoles, that were created during the symmetry breaking, were so diluted during 
this rapid expansion that the monopole density becomes hardly traceable at present. 

So, it is seen that inflationary models could solve the problems of standard Big Bang 
cosmology more or less satisfactorily. But these models also suffer from some problems, the 
most famous one being the ''graceful exit problem" . The problem originates from the fact that 
we see the galaxies and other structures around us. The universe must have a decelerated 
expansion at some epoch of time in order to facilitate the galaxy formation. This demands 
that the universe has to come out of this inflationary phase. This problem that how the 
universe comes out of this rapid expansion phase and enters a decelerated phase of expansion 
is termed as the "graceful exit problem" . A number of models have been suggested for the 
solution of this problem. All of them have their own merits and pitfalls the details of which 
are not discussed over here. 

So all the major problems in Standard Big Bang Cosmology were believed to be related to 
the early phase of the history of the universe. The present state of affairs in the universe were 
presumably competently taken care of, excepting some flner details like that regarding the 
structure formations. The recent observation that the present universe is accelerating came 
as jolt and thus the late time behaviour also warrants serious theoretical attention. 
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1.4 Observational Evidence of the Present 
Acceleration : 

The evidence in support of an accelerating universe stems from the observations of the 
luminosity- redshift relation of type la supernovae. 



Source 



O 

Observer 



Figure 1.2: Measurement of luminosity distance dL from absolute and apparent luminosities 

In a static universe, if one considers a luminous object S emitting a total power L ( also 
called absolute luminosity ), then the intensity I ( called apparent luminosity ) detected by 
an observer at O at a radial distance di from the luminous object (as shown in Figure 11.21) is 
given by 

The quantity 



is known as the luminosity distance. In a static universe, the luminosity distance is equal to 
the actual distance. In an expanding universe however, the intensity detected by the observer 
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gets reduced because the energy of a photon emitted gets redshifted due to the cosmological 
expansion [15]. Because of this expansion, the detected energy gets reduced by a factor of 

-^ = l + z 1.21 
a[t) 

where a{t) is the scale factor of the universe at some cosmic time t, to is the present time and 
z is the redshift parameter , given by 2; = where A is the emitted wavelength and the 
wavelength received is A + AA. 

Thus in an expanding background, the observed apparent luminosity can be written as 

A7ra{tofx{z)\l + zf ^ ' 

where x{z) is the comoving distance of the luminous object, expressed as a function of the 
redshift z. This implies that in an expanding universe, the luminosity distance dL{z) is related 
to the comoving distance x{z) as 

dL{z) =x{z) il + z) (1.23) 

where a is normalized such that a (to) = 1- 

As the light geodesies in a spatially flat expanding universe obey the relation 

cdt = a{z) dx{z) , (1.24) 



one can eliminate x{z) using equation (ll.23p and express the expansion rate of the universe 
if = ^ in terms of diiz) as 

dz \{l+z)) 

Thus if the absolute luminosity of a distant object is known, its apparent luminosity can be 
measured as a function of z and from equation (11.20p . the luminosity distance di can be 
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calculated as a function of redshift z. The expansion history H{z) can then be deduced by 
differentiating equation (]1.25p with respect to z. On the other hand, if a theoretically predicted 
H{z) is given, the corresponding diiz) can be predicted by integrating equation fll.25p as 

Then this predicted di^z) is compared with the measured di^z) to test the consistency of a 
model. In practice, the astronomers do not use the ratio of absolute over apparent luminosity. 
Instead they use the difference between apparent magnitude m and absolute magnitude M 
given by the relation 

m - M = 2.5logio (y) • (1-27) 

It is important to see how the deceleration parameter q, defined by equation (11.91) . can be 
estimated from this observed luminosity at a given redshift. For small values of the time scale 
and the distance compared to the respective Hubble scales, the scale factor a{t) and hence z 
can be written as power series, 

rT2u 4. \2 



z = Hoito - h) + (^1 + f j H^ih - hY + , (1.28) 

where the quantities with a suffix zero indicate their present values. Using this expression, 
the luminosity distance di can be written as 

1 



dL = 



z + -{l-qo)z^ + 



;i.29) 



and hence the apparent luminosity I becomes 



2 



The apparent and absolute magnitudes will then be related as 
m - M = 25 - 5logioHo (km/sec/Mpc) + 5logiocz {km/ sec) + 1.086(1 - qo)z + (1.31) 
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The apparent luminosity is smaller for a negative go and thus the dimmer appearance of the 
supernovae calls for a negative go- The expansion history of the universe is very well depicted 



past 




Accelerating 
Universe 



Luminosity 
Distance 



furtlier 
away 



now 



redshift z 



past 



Figure 1.3: Hubble diagram for an accelerating universe 



by the Hubble diagram where the x-axis shows the redshift z of a luminous object and the 
y-axis shows the physical distance Ar to those objects. This Ar is in fact the luminosity 
distance. As the redshift z is related to the scale factor a(t) at the time of emission of 
radiation via equation fll.2ip . whereas Ar is related to the time in past when the emission was 
made, therefore the Hubble diagram provides information about the time dependence of the 
scale factor a{t). 

The slope of this diagram at a given redshift denotes the inverse of the rate of expansion 



In an accelerating universe, the slope of the Ar vs. z curve is larger at high redshift. 
Thus, at a given redshift, luminous objects appear to be at a greater distance, i.e, dimmer as 
compared to an empty universe expanding with a constant rate (see Figure [L3|l . 



H{z), i.e. 



Ar 



1 





c z 



Observational Evidences 17 

In the construction of Hubble diagram, those luminous objects are used whose absolute 
luminosity is known and therefore by measuring their apparent luminosity, their distances 
can be calculated. Those luminous objects are called standard candles or distance indicators. 
Type la supernovae serve as excellent standard candles for estimating the luminosity distance 
(11- Type la supernovae are explosions believed to occur in binary star systems where one 
of the companions has a mass below the Chandrasekhar limit 1.4Mq and thus become a 
white dwarf supported by degenerate electron pressure after hydrogen and helium contained 
in the star are burnt up. Once the other companion reaches the red giant phase, the white 
dwarf starts accreting matter from the companion star. Once the mass of the white dwarf star 
becomes equal to the Chandrasekhar limit, the gravitational pull overcomes the degeneracy 
pressure and the white dwarf starts to shrink. This increases the temperature and results in 
the carbon fusion. This leads to violent explosions, commonly called supernovae explosions. 
These explosions are detected by a light curve whose luminosity increases rapidly in a time 
scale of less than a month, reaches a maximum and disappears in a timescale of 1-2 months (see 
figure 11.41) . Type la supernovae are characterised by the absence of hydrogen and abundance 
of silicon in the spectrum. 

Supernovae are preferred standard candles mainly for the following reasons : 

i) These objects are highly luminous. This high absolute luminosity of supernovae la (M = 
—19.5) ensures that they can be seen from large distances ( ~ 1000 Mpc ) and thus are useful 
for measuring various cosmological parameters. 

ii) The dispersion in supernovae luminosity at maximum light is extremely small and the 
corresponding change in intensity is ~ 25%. 

iii) Their explosion mechanism is fairly uniform and well understood. 
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Figure 1.4: Light curve for a typical Supernova la 

However, the major problem in using Sn la as standard candles is that they are rare events 
- for instance, in our galaxy they occur only a few times in a millennium. Also it is not easy 
to predict a supernovae explosion. However, the key feature of the high precision tools of 
observation is that one can now detect supernovae in other galaxies also. The basic strategy 
employed in observing supernovae la is as follows [HI [TTl [HI [191 (201 [21] : 

i) A number of wide fields of apparently empty sky are observed. With modern instruments 
on a 4 meter-class telescope, tens or thousands of galaxies are observed in a few patches of 
sky. 

ii) three weeks later, the same galaxies are observed once again. 

iii) The images are then subtracted to observe the supernovae explosions. 

The result of this observation strategy is a set of Sn la light curves in various bands of spectrum. 
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From the light curves, their peak apparent luminosity is used to construct the Hubble diagram. 

The first project in which supernovae were used to determine the energy associated with 
the cosmological constant was carried out by Perlmutter et al. in 1997 [22]. This project was 
named as Supernovae Cosmology Project (SCP). In one year they discovered seven supernovae 
at redshift 0.35 < z < 0.65 and observed them with different telescopes from the earth. 
The Hubble diagrams they constructed were in good agreement with a standard decelerating 
Friedmann cosmology. However, one year later they updated their results by including the 
measurements of a very high redshift {z ~ 0.83) Supernovae la [TB|. This dramatically changed 
the scenario and a decelerating universe was ruled out at about 99% confidence level. This 
result was confirmed independently by another pioneer group, viz, High-z Supernovae Search 
Team (HSST) [T7j. They had discovered 16 supernovae at a redshift 0.16 < z < 0.62 and the 
results indicated an accelerating universe at a 99% confidence level. 

In 2003, Tonry et al. [IH] reported the results of eight newly discovered supernovae in the 
range 0.3 < z < 1.2. Their results reinforced the previous findings of accelerated expansion 
and also gave the confirmation of decelerated expansion at 2; > 0.6. So, obviously the universe 
must have had a transition from deceleration to acceleration in the past. This transition was 
confirmed and pinpointed by Riess et al. in 2004 [13]. They included 16 new high redshift 
supernovae and after analyzing all available data, constructed a reliable and robust data set 
consisting of 157 points which is known as Gold data set. With this new data set, they could 
clearly identify the transition from decelerated to accelerated expansion at z ^ 0.46 ± 0.13. 
However, it was not easy to conclude whether the data favour an accelerating or decelerating 
universe from only the Hubble diagram corresponding to the Gold data set because of the error 
bars present in the data set. This would be easier if the Hubble diagram of figure [T73] where the 
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distance is plotted vs redshift is superposed with the distance- redshift relation d^j^^^^{z) of an 
empty universe with H{z) constant. So, a more efficient plot was used for this purpose, viz, the 
logarithmic plot of diiz) vs. <f^'^^^{z) which could easily distinguish between an accelerated 
and decelerated expansion. Such a plot for the Gold data set clearly indicated that the best 
fit is obtained by an expansion which was decelerated at the earlier times {z > 0.5) and 
accelerated at recent times {z < 0.5). Attempts have been made to explain the observed 
dimming of supernovae at high redshift by considering that this apparent dimming is due to 
the scattering of light by intergalactic dust or grey dust or even due to evolution of Sn la. 
However none of them has a &m footing [T3| [2^ [25] and thus strengthens the belief that the 
universe at present is undergoing an accelerated phase of expansion. This is also conffimed 
by the highly accurate Wilkinson Microwave Anisotropy Probe (WMAP) data [3 El El El E] • 



1.5 Search for the Dark Energy : 

This observed acceleration of the universe brings in trouble as gravity is attractive and this 
acceleration can not be driven by the attractive gravitational properties of regular matter. 
So, obviously an additional component is required which can give rise to an effective repulsive 
gravity so that matter can move away from each other with an acceleration. The obvious 
question to address is therefore, ''What should be the properties of this additional component^" 
The answer can be obtained by comparing Newtonian gravity with Einstein's gravity. 

In Newtonian gravity, the acceleration of a test particle having mass m under the influence 
of gravity is given by 

md = 
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which gives 



a 47rG' 

- = —p (1.33) 



as M = |7ra'^p. 



On the other hand, in Einstein's gravity, from equations (II .41) and (ll.Sp one can obtain 

d AttG , . , , 

- = ^(P + 3p). (1.34) 



This can be written in terms of equation of state parameter w as 

- = —Pil + Sw) (1.35) 

a 3 

which gives a measure of the acceleration of the universe. 

So, from equations (11.33^ and ( ]1.34p . it is evident that unhke that in Newtonian gravity, in 
Einstein gravity, pressure also plays a major role along with the density in determining the 
space-time dynamics of the universe. In other words, "pressure carries weight in Einstein's 
gravity". 

Now, as the universe at present is accelerating, a > 0. From equation (11.341) . a positive d can 
be obtained only if p is sufficiently negative such that p + 3p < or equivalently w < — |. So 
in order to explain the observed acceleration of the universe there is indeed the requirement 
of some form of matter which can generate sufficient negative pressure. This particular form 
of matter, now popularly referred to as ^^dark energy" , is believed to account for as much as 
70% of the present energy of the universe. 

A large number of possible candidates suitable as dark energy component have appeared 
in the literature. None of them has a clear advantage over all others. All have their merits, 
but none perhaps has a firm theoretical footing. Its effect has been, as discussed, is to provide 
a sufficient effective negative pressure. Nothing is known about its distribution vis-a-vis the 
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dark matter, except that it does not cluster at any scale lower than the size of the universe. 
A few of these candidates are described below. 



1.5.1 Cosmological Constant Models :- 

The simplest dark energy candidate is the cosmological constant A introduced by Einstein in 
1917. With the introduction of the A-term, the Einstein's equation gets modified as 

G^, = 87cGT^, + Ag^, . (1.36) 

Einstein originally introduced the A-term on the left hand side of the field equation in order to 
obtain a static universe. But later on, when Bubble's observations suggested that the universe 
is expanding, he himself rejected the A-term. However, A was again brought into being in 
early 1980's with the inflationary model of the universe p3|. During inflation, the universe 
was dominated by the vacuum energy p^ac as discussed earher. As the equation of state for 
such energy is 

Pvac Pvac 0, 

both p^ac and p^ac are constants, and the Einstein field equations will effectively look like 

Qj Jv , , , 

3^ + 3^ = A, 1.37 
2- + ^^ = A, 1.38 

where p^ac = -Pvac = g^- 

Here A is not introduced arbitrarily, but rather attains the significance of the vacuum energy 
density. By solving these equations one gets an exponentially expanding model of the universe 
at a very early epoch which provides solutions to the problems of Big Bang theory. As the 
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simple inflationary model had its own problems, such as that of the 'graceful exit', or the huge 
discrepancy in the theoretically predicted value of A and the one suggested by observations, 
other more complicated models took over and A had to hide itself into oblivion. In the late 
90 's, with the observation that universe is at present accelerating, A again came back strongly 
after a brief period of hibernation. 

For an FRW universe, the Einstein's equations with a cosmological constant take the form 

3— + 3— = SvrGp + A, (1.39) 

d a"^ + k ^ . / N 

2- + —— = -SnGp + A. 1.40 
a 

From equations (11.391) and (11.401) . one easily arrives at the relation 

a AnG , ^ , A /, 

- = ^(P + 3p) + - , (1.41) 

a 6 6 

from which it is clear that A is capable of providing an acceleration (a > 0) if 

A > 47rG(p + 3p) . 

In a very simple case, when the universe is spatially fiat {k = 0) and is dust dominated (p = 0), 
equations (11.391) and (11.401) gives an exact analytic expression for the scale factor as 

2/3 

2V3*' 



a{t) oc I sinh-W-t 1 . (1.42) 



For a very small t (early epoch), a{t) oc t'^^'^ and hence q = which gives a decelerated 
expansion in the early phase of dust dominated era as expected [26]. On the other hand, 
for large t (i.e, present epoch), a oc eV^* which gives an accelerated expansion. This model 
is popularly known as ACDM model. However, there is a major problem related to the 
cosmological constant A. Field theory suggests that the lower limit of the value of vacuum 
energy density should be 
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p,ac ~ 10^ GeV^ , 

whereas the current observational upper hmit on the cosmological constant is 

Pvaclo ~ 10-2^ g/cm^ ^ 10"^^ GeV^ , 

where the subscript '0' stands for present time. So, there is a discrepancy between the predicted 
value and the observed value by 10^^ orders of magnitude. The cosmological constant A is 
expected to have a large value during early epoch so as to resolve - via inflation - the horizon 
and flatness problems; at the same time it requires a low value at the present epoch to avoid 
conflict with the observations. This problem is referred to as the ^^cosmological constant 
problem" because of which A has lost some of its ground. 

To avoid this problem, some dynamical models of dark energy have been proposed where 
A, instead of being a constant, is a slowly varying function of time. 

1.5.2 Varying A Models :- 

The first candidate for a dynamical model of dark energy was the time varying cosmological 
parameter A[t). A number of phenomenological models have been described to introduce this 
dynamical A-term as the new form of matter [271 1211 [291 EOl EH |32] and is taken as a function 
of the scale factor a{t) or the cosmic time t. However, there is no strong physical motivation 
behind this choice. In fact, the dynamical A-term leads to some problems. If A is considered 
as a constant, then the introduction of the term Ag^^i^ in the Einstein field equations does 
not affect the conservation equation. It is so because the conservation equations come as a 
consequence of Bianchi identity 

G(t- = , 
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which remains the same as for a constant A, 



iAgn;u = A,.g^'' + Agtl^ = 0. 



However, in a varying A model, the Bianchi identity leads to a modified conservation equation. 
If the normal matter is allowed to satisfy its own conservation equation of the form given in 
equation (11.61) . A automatically turns out to be a constant. The way out is to consider some 
interaction between A{t) and matter such that one grows at the expense of the other. But in 
order to incorporate this interaction, one has to sacrifice the matter conservation equation. In 
most of the cases, the choice of the form of interaction does not have any physical background. 
So, the second natural choice for a dynamical A-term is a scalar field having some potential 
V{(f)) and is popularly called the ^^quintessence scalar field" . 

1.5.3 Quintessence Models :- 

Inspired by the role of a scalar field with a suitable potential in providing an infiationary 
scenario in the early universe, attempts are made to formulate a "dynamical dark energy" 
model based on scalar field cosmology. With the introduction of a scalar field minimally 
coupled to gravity, the action gets modified as 





For a spatially fiat FRW cosmology, Einstein's 



field equations take the form 




+ -0' + V{(j)) , 

= -Pm - ^0' + V{<P) 
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Variation of the action with respect to the scalar field leads to the wave equation for as 



a 



(j) + 3-(f) + V'{(j)) =0 (1.46) 
ct 

where a prime indicates differentiation with respect to 0. 

In view of the equations fl 1.441) - fll.46p . the matter conservation equation 

Pm + 3-{pm+Pm) = (1.47) 

a 

is not an independent equation, and comes as a consequence of the Bianchi identities. Thus 
one is left with five unknowns, a, pm, Pm, 4> and V{(j)) and only three independent equations 
to solve them. If the equation of state Pm = Pm{Pm) and the form of potential V = V{(f)) 
are given, the system of equations is closed. Equations fll.44p and fll.45p reveal that the 
contribution to the effective energy density and pressure from the scalar field is 

P<p = \4>' + V{ct>) , (1.48) 

P4> = - V{^) , (1.49) 

where the first term on the right hand side represents the kinetic part whereas the second 
term represents the potential part. 

In the matter dominated era the fiuid pressure Pm = 0, and equation (11.341) yields the result 



a 
a 



AttG 
3~ 

AtcG 



[pm + Pcj, + 3p</,] 



Pm + 



2Vi<P) 



;i.5o) 



If 1^(0) grows to a sufficiently large value such that 



2V^(0) >p^ + 202 , 



(1.51) 
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d/a has a positive value consistent with the present observation. These models have a further 
interesting possibility that if V{(f)) has a smaller value compared to ^""^^"^ in the early matter 
dominated era and grows later to dominate the dynamics, the model would exhibit a signature 
flip in q from a positive to a negative value at a recent past. A fine tuning of the parameters of 
the theory and the arbitrary constants of integration might lead to models consistent with the 
details of observational data. These features have made scalar field cosmologies a very active 
area of research, and now the scalar field with a potential giving rise to a negative pressure 
is specifically called a '^quintessence field" and the generic term for any field, that drives the 
late surge of accelerated expansion of the universe, has become 'dark energy\ 

The system of equations suggest that if a particular dynamics of the universe, i.e, the 
temporal behaviour of the scale factor a{t) is chosen, it is possible to find a potential consistent 
with that. Amongst the host of quintessence potentials available in the literature, some give an 
ever accelerating model, some give an acceleration in the large t limit, and a few indeed allows 
for a smooth transition from a decelerated to an accelerated expansion in the matter dominated 
era itself consistent with the observation that the universe has entered the accelerated phase 
of expansion in the late matter dominated era. 

One typical example of the latter kind of a quintessence field is the one introduced by Sen 
and Sethi [33], where the potential is a double exponential, given by 

= «(e2^'^ + e-^^-^) + V^o • (1.52) 

where a, [3 and Vq are constants. The solution of the scale factor comes as sin h'^(t/to) . For 
a small t, the solution for the scale factor is like ~ whereas for a high value of t, a ~ e^*/*°. 
For /3 < 1, the model has a decelerated expansion in the early epoch, but indeed accelerates 
at a matured age. Similar kind of potentials have been used by Barreiro et al [33] and Rubano 
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et al [35] as well. 



A particularly interesting class of quintessence fields are the so called ^Hracker fields" for 
which the potential is steep enough to satisfy the condition 



All potentials satisfying this condition lead to a common evolutionary path for the fields 
starting from a wide range of initial conditions [36] • The energy density of the tracker fields 
has an evolution which mimics that of the background matter during most of its history 
(see Figure ll.Sp and dominates over the matter density at later stages of the evolution to 
generate acceleration. Thus, these models can alleviate the coincidence problem which poses 
the question why the dark energy sector should dominate the dynamics of the universe at 
present. 






z+1 



Figure 1.5: The quintessence field rolling down an inverse power law potential tracks first 
radiation, then matter and finally dominates the energy density of the universe at present. 
(From Zlatev, Wang and Steinhardt JE^) 
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Some simple examples are the inverse power law potential 



yi<p) = (1-54) 



with some restriction on a, and the exponential potential [37l [38l [39] 

V{<j)) = Vo e-^'^ . (1.55) 

Within the tracker framework, a potential with a versatile ambition was given by Sahni and 
Wang [lO]. This potential is given as 



V{(P) = K)[cosh(A0) - 1]^ . (1.56) 
For I A0 I >> 1 and < 0, the potential is exponential, 

V{(f)) oc exp(— pA0) , 

and W0 = ^ is very close to Wm = 

On the other hand, when | A0 | << l,V{(f)) oc (A0)^^. Then, the average equation of state 
parameter is given by 

< >= 

For p < |, the model describes quintessence whereas for p = 1 it can play the role of a cold 
dark matter. So, this model is able to describe both dark matter and dark energy within a 
tracker framework (see [ITl H2] ). 

Another useful potential with interesting features was proposed by Zlatev et al. ^36j as 

V{(j)) = Vo \e^^^/^ - l] . (1.57) 
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The advantage of this potential is that it can significantly alleviate the fine tuning problem 
and P(j) can come to dominate over the present matter density from a large number of initial 
conditions. 

There are of course a lot of other examples. However, despite the many attractive features 
of these quintessence potentials, a degree of fine tuning does remain in fixing the parameters 
of the potentials. The generic problem of quintessence models is that the potentials are all 
taken arbitrarily and none of them has a sound physical basis. But indeed if we know the kind 
of acceleration required, more often than not, one can find out the form of potential which 
generates the desired acceleration. 

1.5.4 Non-minimally Coupled Scalar Field Models :- 

In these models, the scalar field is non-minimally coupled to gravity such that the general 
action is of the form 



where and g{(f)) are some functions of the scalar field 0, Go is the Newtonian constant of 
gravity. 

Comparing this action with the Einstein's action, one gets 



So, in this type of scalar field models, G is not a constant but is some function of the scalar 
field 0. 

The simplest of the non-minimally coupled scalar field models is the Brans-Dicke theory. 
This theory was first introduced to incorporate Mach's principle in a relativistic theory of 
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gravity. According to this principle, the inertia of a body is not an intrinsic property of its 
own, it rather depends on the mass distribution of the rest of the universe. 



In Brans-Dicke theory 



= 



9i<P) = T 



1.59) 



uj being a dimensionless constant parameter. With these, the action for the Brans-Dicke theory 



takes the form 



5 = y ^ 



(j)R UJ 



;i.6o) 



IGtvGq (j) 

If we consider the matter field to consist of a perfect fluid, then the field equations for a 
spatially flat Robertson - Walker spacetime become. 



a" 



a 

2- + - 

a 



2 ^ Pm 

-SttGo— - --^ - 2-- 



;i.6i) 
;i.62) 



Also, the wave equation for the Brans-Dicke field is 







a ■ SttG, 
+ 3-6 = ^ 



1.63) 



Dicke [H] in 1962 framed an alternative version of their scalar tensor theory [13] by a simple 
redefinition of units. By effecting a conformal transformation 



;i.64) 



the action given by equation fll.60p becomes. 



S 



R 



16nG 



(2cu + 3) ^ ^ „ - 



1.65) 
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where an overhead bar represents quantities in new frame and \1/ = ln0. 

Comparing this action with equation (]1.60p . it is seen that in this version (f) is no longer coupled 

to R. So, the effective constant of gravitation G, which is a function of the scalar field as 

G = ^, (1.66) 

in the original version of the theory, now becomes a constant. 

The field equations f ll.Gip and (11.621) in the new frame look like 



= p„. + (1.67) 



a a (2uj + 3) • o , , 

2^ + - = -Pm- ^ , , 1.68 
a 4 

and the wave equation becomes 

^ a^^Anj^^ (1.69) 

a (2cu + 3) ' ^ ' 

where an overhead bar indicates quantities in new frame and = (pa^. The equations are 
written in units where SttGq = 1. The density and pressure of the normal matter in this 
version are related to those in the original version as 

Pm = (t)~'^Pm, and pjn = (t)~'^Pm- 

The resulting field equations in the new frame look more tractable than the original version. 
Furthermore, the equations in this version give insight regarding the comparison of the energies 
of different components of matter. For example, it is clearly seen from equation (ll.67p that 
the contribution to the energy density by the scalar field is given by 
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However, one has to pay some price for it. In the transformed version although G becomes a 
constant, the rest mass of a test particle becomes a function of the scalar field [44] and one 
has to sacrifice the equivalence principle. So, the geodesic equations are no longer valid and 
indeed the physical significance of different quantities in this version of the theory is somewhat 
obscure. However, because of its computational simplicity, it is easier to arrive at some 
solutions in the transformed version. And the problem can be resolved by transforming back 
to the original atomic units where one can talk about the various features more confidently. 

Although General Relativity (GR) is a better theory of gravity and it enjoys experimental 
evidences in its support, for various reasons Brans-Dicke (BD) theory continues to enjoy an 
alive interest. One important advantage of BD theory is that it becomes indistinguishable 
from GR in the limit uj ^ oo. This limit has now shown to have only a restricted application 
[15] . but the PPN parameters calculated in BD theory [16] clearly shows that at least in the 
weak field limit the predictions from local astronomical observations in this theory will be 
same as that in GR in the large u limit. For these reasons, there is a popular notion that 
BD theory is perhaps the most natural generalization of GR. However, the local astronomical 
observations suggest that if BD theory has to be consistent, then u ~ 10'^ which renders the 
theory practically indistinguishable from GR 

BD theory proved to be useful in providing clues to the solutions for some of the outstanding 
problems in cosmology. In 1981, Guth [H] proposed the inflationary model of cosmology. This 
model could provide solution to many of the cosmological problems as discussed earlier, but it 
suffered from the ^graceful exit prohlem\ A large number of models were introduced to solve 
this problem having their own merits and demerits. 
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In 1984, Mathiazhagan and Johri [18] addressed the problem in Brans-Dicke (BD) theory 
[43] . Under this framework, it was shown that along with a vacuum energy, the scale factor 
grows as a power function of time. Using a similar technique. La and Steinhardt [19] presented 
the extended inflation model" in order to get a sufficient slow roll of the scalar field so 
that there is sufficient time for the completion of phase transition and thus the graceful exit 
problem could be resolved. However, this leads to unacceptable distortions of the microwave 
background |.50j. To solve this problem, Steinhardt and Accetta [5T] developed the ''hyper- 
extended inflation model" . Later Brans-Dicke theory was used for finding a solution to the 
graceful exit problem with a large number of potentials [52] , where the infiaton field oscillates 
during the later stages of evolution and the universe comes out of the inflationary phase. 

BD theory has also found applications in solving a few recent cosmological problems, such 
as, the quintessence problem. A number of models have been presented where the Brans- 
Dicke scalar tensor theory could potentially solve the problem of quintessence as it leads to 
non-decelerating solutions for the scale factor in the present matter dominated universe. In 
some of these models [53l the BD theory is modified by incorporating a potential V{(f)) 
which is a function of the BD scalar field itself which could drive the acceleration. However, the 
problem with these models is that V{(f)) is put in by hand and there is no physical motivation 
behind the choice of form of V{(j)). 

A few other models have also appeared in the literature where the cosmic acceleration is 
obtained with a quintessence field in BD theory [55j. However, this result hardly provides 
any improvement on the corresponding GR result. Recently Banerjee and Pavon ^56j have 
proposed a model where the BD theory could explain the present accelerated expansion of the 
universe without resorting to a cosmological constant or quintessence matter. This is better 
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in the sense that one does not have to invoke any additional quintessence field to explain the 
acceleration. However, this model have problems in providing a decelerated expansion in the 
radiation-dominated epoch. 

The general defect of all these models is that none of them could provide a smooth transition 
from the decelerated to accelerated phase of expansion and they rather provide an acceleration 
in some limit. Also in all of these models, a consistent accelerated solution is obtained only 
for small negative values of the Brans-Dicke parameter u. This is in sharp contrast to the 
value obtained from local astronomical experiments which predict the value of u to be of the 
order of a thousand [57] . Attempts have been made to overcome this problem by considering 
a modified version of Brans-Dicke theory, called Nordtvedt's theory, where the parameter uj 
is a function of the Brans-Dicke scalar field instead of being a constant [57J. In this case the 
field equations fll.6ip and (11.621) remain intact, but the wave equation fll.63p gets modified as 



Bartolo and Pietroni [SH] have pointed out that a varying uj can indeed explain the late time 
behaviour of the universe. Also Banerjee and Pavon [56] showed that a varying u theory could 
give rise to a decelerating radiation model followed by an accelerating model in the matter 
dominated universe. However, none of these could provide smooth transition from deceleration 
to acceleration in the matter dominated era itself. So, a thorough survey of varying u theory 
is indeed warranted to check if it gives rise to a model of the universe which can explain the 
transition from decelerated to accelerated phase of expansion in the matter dominated epoch 
itself with some high values of uj consistent with local astronomical experiments. 
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1.5.5 Curvature Driven Accelerating Models :- 

The scalar field models or the cosmological constant models are amongst the most popular 
candidates of dark energy component. Recently, an attempt along a different direction is also 
gaining attention. This effort explores the possibility of whether geometry by itself can serve 
the purpose of providing late time acceleration of the universe. 

The idea actually originates from the experience of inflationary models. It was shown 
by Starobinsky [59] and Kerner et al pOl [M] that higher order modifications of the Ricci 
curvature R, in the form of or R^^R^" in the Einstein - Hilbert action, could generate 
sufficient acceleration in the very early universe. However, with the evolution of the universe, 
R is expected to fall off. This leads to the question whether the inverse powers of i?, which 
becomes dominant during the late time, can help driving the recent acceleration. 
The action gets modified as, 

1 



S 



-J{R) + i« 



^(fx (1.71) 



.167rG'" 

where the usual Einstein - Hilbert action is generalized by replacing R with an arbitrary 
function f{R). A variation of this action with respect to the metric yields the field equations 



as 

1 

2 



G^u R^j,u c^Rq^iu ^^1/ + , (1.72) 



where T^^ represents the contribution from the curvature and denotes the energy momen- 
tum tensor components for the matter field scaled by a factor of j^^- Here the choice of units 
SttG = 1 has been made. 
T^^ is explicitly given as, 

1 n 1 

(1.73) 



^^^^ f'{R) 



^9tiu{f{R) - Rf'{R)) + fiRy^^ig^aQup - gf^uQap] 
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where a prime indicates differentiation witli respect to tlie Ricci scalar R. 



For a spatially fiat Robertson- Walker spacetime, where 



ds'^ = dt^ - a^{t) dr"^ + r'^dO'^ + r^sin^ Odcj)"^ 



(1.74) 



the field equations (11.721) take the form 




r{R) 1 2 




[/(i?)-i?/(i?)]-3-i?r(i?)|+p. 



m 



(1.75) 




1 



{ 



2-^Rf\R) + Rf\R) + R^r{R)-]^{f{R)-Rf\R))^-p^ . (1.76) 



It is evident that if f{R) = R, the field equations fll.75p and fll.76p take the form of usual 

Einstein field equations. 

The Ricci scalar R is given by 



one actually has a system of fourth order differential equations. Depending on the functional 
form of f{R), some of the terms on the right hand side of equation fll.76p can provide an 
effective negative pressure and generate sufficient acceleration. 

A substantial amount of work has already been done along this line by choosing various 
functional forms of f{R)- Capozziello et. al [621 [63] considered f{R) = R"' and showed that 
it leads to an accelerated expansion for n = —1 and n = |. The dynamical behaviour of i?" 
gravity has been studied in detail by Carloni et. al [7j. Carroll et. al [5] used a combination of 
R and in the action and a conformally transformed version of the theory where the effect of 
curvature is formally taken care of by a scalar field having some potential. They showed that 
it could generate a negative value for the deceleration parameter q. Vollick [6j used term 



R = -6 - + ^ 

a 



(1.77) 



which involves a second order derivative of the scale factor a. As equation fll.76p contains i?, 
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in the action and obtained an exponentially expanding and hence accelerating model for the 
universe. It deserves mention that Vollick actually employed a Palatini variation, so the field 
equations are different from equations fll.75p and (11.761) . Nojiri and Odinstov [8j considered 
the Lagrangian of the form 



and showed that it is indeed possible to obtain an infiation at the early stage and a late time 
accelerated expansion from the same set of field equations. Other interesting investigations 



acceleration. However, all these models mentioned above have problems regarding the stability 
Furthermore, most of them either resort to a piecewise solution for large R and small R, 
or provide acceleration in some limit or an eternally accelerating model. But none of them 
could show the transition from decelerated to accelerated phase of expansion in the same 
matter dominated regime. But still these investigations open up an interesting possibility for 
the search of dark energy in the non-linear contributions of the scalar curvature. 

1.5.6 Chaplygin Gas Models :- 

A chaplygin gas model is also one of the important candidates for solving the dark energy 
problem. The Born-Infeld lagrangian density 



L = R + R^ + i? " where m, n are positive integers. 



include the choice of f{R) as sinh ^{R) [68] or Ini? [lOj, which also could provide late time 





leads to the chaplygin gas obeying the equation of state 



A 



(1.79) 
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where A = Vq^ . 

The chaplygin gas can also be derived from a quintessence lagrangian 

L = \^'- V{ct>) 

with the potential [71] 

2 \ cosh 



y(0) = ^ cosh30 + — . (1.80) 



The conservation equation 



dipa'') + p dia-") = (1.81) 



immediately gives 



P=\A + -^ (1.82) 



a6 



where 5 is a constant of integration if the equation of state is given by equation fll.79p . 
The more general form of equation fll.79p leads to the equation of state for the generalized 
chaplygin gas given by (see [72] and references therein ) 

p=-— (1.83) 

where < a < 1 and A is a positive constant, a = 1 gives back the old chaplygin gas model. 
In the framework of FRW cosmology, this equation of state yields solution of the Einstein 
equations and leads to density evolving as 

where a is the scale factor of the universe and B is an integration constant. From equations 
(11.821) and fll.84p . it is clear that for small t, i.e, for small a, one can obtain a dust dominated 
model (^p ~ whereas for large t, p \fA or (p ~ and it behaves like a cosmological 

constant and provides acceleration. 
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1.5.7 Phantom Dark Energy Models :- 

Caldwell [73] pointed out that a very good fit to the luminosity-distance curve (Figure ll.3p 
can be provided by a dark energy component which violates the weak energy condition so 
that the equation of state w < —1. He dubbed this candidate as ''phantom dark energy" . A 
study of high-z supernovae [21j also reveals that the dark energy equation of state has 99% 
probability of having a value < —1 if no constraints are imposed on Qm- 

A number of models appeared in the literature where the dynamical nature of phantom 
energy was constructed by taking a kinetic term with a 'wrong' sign in equation fll.49p so 
that it can give rise to the present acceleration of the universe [TH [751 [76l [77]. Phantom dark 
energy models have also been studied in Brans-Dicke theory [78] or in an interacting scenario 
where the phantom field is coupled to some other field [791 [HH] • However, these models suffer 
from the problem of instability at the quantum level ^1] as there is no proper ground state 
because of the negative kinetic energy. Also models with w; < — 1 suggest that the effective 



velocity of sound in the medium v = y \ dp/ dp \ can become larger than the velocity of light. 
The phantom models imply a pathological behaviour for the cosmological model at a finite 
future. If teq denotes the time when matter density and phantom energy density become equal, 
then the scale factor of the universe grows as 



a{t) ^ a{teq) 



^ 2 

[l + w) w 

teq 



3(l+w) 

, w <-l . (1.85) 



Therefore, when t (^Ih^) ^eq, (i{t) — »■ c>o , i.e, the scale factor diverges in a finite time. 
At that epoch, Hubble parameter H also diverges implying that the expansion rate of the 
universe reaches an infinite value in a finite time. This situation is termed as 'Big Rip\ Thus 
the universe dominated by phantom energy culminates to a future curvature singularity ( 
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See [HSllHSllTSllHaiHnilHHllHTlEHlEnil^ However, some other models have also been 

investigated where w < — 1 is attained without considering a negative sign for the kinetic 
term. These models are called 'Braneworld models' p2l [15] which has Wejf < —1 today, but 
does not run into a 'Big Rip' in finite future. 

1.5.8 Braneworld Models :- 

Braneworld cosmology suggests that we could be living on a four dimensional '6rane' which is 
embedded in a five or higher dimensional ^bulk\ It is considered that matter fields are confined 
to the brane whereas gravity is free to propagate throughout the bulk (for a comprehensive 
discussion, we refer to the lectures by Roy Maartens [M]). In the Randall-Sundrum (RS) [95] 
scenario, the equation of motion of a scalar field propagating in the brane is given by 

^ + 3H<j) + V'{(f)) = , (1.86) 

where 

P = l<P' + Vi^) ■ (1-88) 
Here, e is an arbitrary constant and cr is the brane tension which relates the four-dimensional 
Planck mass (m) and the five- dimensional Planck mass (M) as 

Also the four-dimensional cosmological constant A4 on the brane and the five-dimensional 
cosmological constant Af, on the bulk are related as 
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Equation (11.871) contains an additional term ^ because of which the damping experienced 
by the scalar field as it rolls down the potential dramatically increases so that infiation can 
be sourced by potentials, such as V oc e~^^, V oc etc, which are normally too steep to 
produce slow-roll. This gives rise to the possibility that both infiation and quintessence may 
be obtained from the same scalar field. These models are called ^quintessential inflationary 
models^ (see [961 EH EH IHll I100[ llOlj and references therein). An example of quintessential 
infiation is shown in Figure (II. 6p . 



0.8 



0.6 

a 

































Kinetic 1 


Radiative 




; ; Matter : 
; ;Dominated '; 






I 







5 10 15 20 

Logio(a/a^nd) 

Figure 1.6: The post-inflationary density parameter Q is plotted for the scalar field (solid line), 
radiation (dashed line), and cold dark matter (dotted line) in the quintessential inflationary 
model. (From Sahni, Sami and Souradeep [9^]) 

A different way of obtaining an accelerating universe was suggested in the braneworld 
model developed by Deffayet, Dvali and Gabadadze (DDG) |102l 1103] . Here both and a 
were set to zero, while a curvature term was introduced in the brane action so that it takes 
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the form 

S = M'^ R+ R+ L 

matter • 

(1.91) 

J bulk J brane J brane 

The resulting Hubble parameter in the DDG braneworld model is 



2 

where /c = ^ is a new length scale. An important property of this model is that the 
acceleration of the universe is not obtained from any ^dark energy^ component. Since gravity 
becomes five dimensional on length scales R > Ic = 2Hq^(1 — Vlm)^^, it is seen that the 
expansion of the universe is modified during late times instead of early times as in the RS 
model. 

A more general class of braneworld models, which includes RS cosmology and DDG brane 
is described by the action [104[ 1105] 

S = M' [ {R- 2A,) + / {m'R - 2a) + [ L^^tter ■ (1-93) 

J bulk J brane J brane 

For 0" = Afc = 0, it gives back the DDG model, whereas for m = it reduces to RS model. 
Sahni and Shtanov [92j have shown that the Hubble parameter for this action comes out as 



ilmil + zY + + 2^1 T 2J^iJn^{l + zf + Vt, + VLi + , (1.94) 



where ^li = j^jp, = -^^w^ = 3^2^:2 , ^a,, = -g^- (The =F sign refers to two different 



ways in which the brane can be embedded in the bulk (for details see [92]). 

An important feature of the braneworld model given by equation fll.94p is that it can lead 
to an effective equation of state of dark energy We// ^ —1- Also it has been shown that in 
this model, the acceleration of the universe can be a transient phenomenon which ends once 
the universe returns to matter dominated expansion after the current accelerated phase of 
expansion and hence does not fall into the problem of 'Big Rip'. 
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1.6 Outline of the Present Thesis : 



Whatever we have discussed so far indicates that our universe at present is undergoing an 
accelerated expansion preceded by a decelerated one. This means that the deceleration pa- 
rameter q must have a signature flip from a positive value to a negative one in the recent 
past during the matter dominated era itself. This indeed leads to the search for a candidate 
which can drive this transition from deceleration to acceleration. This thesis is a collection of 
papers which investigates some cosmological solutions where this transition is obtained with 
or without the need of "dark energy" . 

The investigations carried out will be presented in four chapters following the introduction. 
The first of them (chapter-2) consists of one paper entitled "Acceleration of the universe with 
a simple trigonometric potential" . In this paper we consider a minimally coupled scalar field 
(j), with some potential V{(j)), as the driver of the late time acceleration. The form of deceler- 
ation parameter q as a function of scale factor a has been chosen such that it has the desired 
property of signature flip and then from this the scalar fleld and required potential are found 
out. 

The flat Robertson- Walker line element is written as 



ds'^ = dt^ - a{tf[dr'^ + r^d^f] . 



(1.95) 



The expression for the deceleration parameter q is taken as, 



Outline of the Thesis 45 

where p is a constant. It is found that for —2 < p < —1, the deceleration parameter q indeed 
has the property of a signature flip. Equation ( 11.96^ integrates to yield 

H = - = A{l + aP), A>0, (1.97) 
a 

where A is an arbitrary constant of integration. 

In the present work, the problem is completely worked out for p = — |. The reason for 
choosing this particular value of p is that it yields q = 0.5 for a very low value of a. For a flat 
matter dominated FRW model without any quintessence matter indeed q = 0.5. So the matter 
dominated universe has a chance to evolve from the standard setting where the formation of 
structures is facilitated. Furthermore, the transition from the radiation dominated era to the 
matter dominated one is quite well understood for q = 0.5 on the matter-dominated side. 

As we are interested in the matter dominated epoch, the fluid is taken in the form of 
pressureless dust and the Einstein equations are written with a minimally coupled scalar field 
having a potential V{(f)). 

With the form of the deceleration parameter assumed in equation fll.96p . the Einstein field 
equations are solved and it is found that the potential V{(f)) is obtained as a trigonometric 
function of the scalar field 0, viz, V{(f)) = ^ cot^ ("^) + 3^^- Also the expressions for 
the dimensionless density parameter Q^n for the normal matter and for the quintessence 
matter are worked out for this model and it is found that the values obtained are well within 
the constraint range obtained from observations [261 El 1107] . 

The equation of state parameters w for the total matter (wt) and that for the quintessence 
matter (w^) are calculated and are found to be consistent with observational requirements. 

This investigation shows that we can have a single analytical expression for q which grace- 
fully transits from deceleration to acceleration and thus has a better footing than the piecewise 
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solutions. Also the form of the potential obtained here is a simple trigonometric function of 
(f) and thus adds to the list of quintessence potentials that serve the purpose of providing a 
presently accelerating universe pUS]. 



The chapter 3 consists of one paper entitled ''Spintessence : A possible candidate as a 
driver of the late time cosmic acceleration" . In this paper we work with a complex scalar field 
model, namely spintessence, proposed by Boyle et al. [109] . The scalar field is taken as, 

^ = 01+^02, (1.98) 

which can also be written as 

^ = 0e^'^* (1.99) 

such that the complex part is taken care of by a phase term. This kind of complex scalar fields 
are already known in the literature with reference to the ^cosmic strings^ [llOj . However, 
cosmic strings have a very specific form of potential whereas in quintessence models, the form 
of potential has to be found out. 

In the present work it has been shown completely analytically that the deceleration pa- 
rameter q indeed has a signature flip from positive to negative values which indicates an early 
deceleration and a present acceleration. In the spintessence model proposed by Boyle et al. 
[109j . it has been considered that u; is a slowly varying function of time which indeed can be 
approximated to be a constant over the entire period of dust dominated era. In the present 
work, we also choose a; to be a constant. This choice along with the field equations leads to 
an equation 

a^^ = -l-mH\ (1.100) 
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where a is the scale factor of the universe, if = ^ is the Hubble parameter and /, m are 
positive constants, which involve the parameters of the theory and constants of integration. 
Now, we make a transformation of time co-ordinate as, 

at a-^ 

which being positive definite indicates that x is a monotonically increasing function of t. We 
now express the various parameters of the model as functions of new cosmic time x without 
any loss of generality or distortion of events. 

With these choices, the Hubble parameter H and the deceleration parameter q comes out 
in terms of the new variable x as, 

H = nt&n{/3 - mnx) (1.102) 

and 

g = -l + ^(^ + m) . (1.103) 
It is evident from equation (11.1031) that q has a zero when 

ai' = ^+m, (1.104) 

the suffix '1' indicates the values of the quantities for q = 0. 

Also from equation fll.1031) it has been shown that ^| ^ is negative definite which indicates 
that g is a decreasing function of a, at least when g = 0. Thus, q definitely enters a negative 
value regime from a positive value at a = ai and H = Hi. So, it has been shown completely 
analytically that a spintessence model can very well serve the purpose of providing the sig- 
nature fiip in q. The important result achieved here is that if q changes its sign in course 
of evolution, this change will be in the right direction, i.e, from positive to negative. Also it 
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has been shown that the flip occurs irrespective of the particular form of potential which is 
a bonus as the form of the quintessence potential is yet to be specified. Another important 
feature of this work is that the value of the scale factor a, where this flip in q takes place, can 
be expressed in terms of parameters of the theory and various constants of integration. The 
latter are arbitrary constants and thus can be adjusted to fit into the observational results. 

The chapter 4 consists of two papers where attempts have been made to explain the late 
time acceleration of the universe in the framework of scalar tensor theories. In the first paper 
we work out the problem in Brans-Dicke theory [13] whereas in the second paper we work in 
Nordtvedt's theory [57] which is a generalized version of Brans-Dicke theory. 

The first paper in this chapter is entitled "^4 late time acceleration of the universe with 
two scalar fields : many possibilities''^ and addresses the ^graceful entry^ problem of how the 
universe transits from a decelerating phase of expansion to an accelerated one. The basic 
motivation for using two scalar fields stems from existing literature on inflation. Mazenko, 
Wald and Unruh [111] showed that a classical slow roll is invalid if the scalar field driving 
inflation is self interacting. Also it has been shown that a single scalar field with a slow 
roll puts generic restrictions on the potentials driving inflation [112] . These problems led to 
the belief that a successful inflationary model requires two scalar fields [113j . The present 
work uses this idea of introducing two scalar fields. One of these fields is responsible for 
the present acceleration of the universe, called the quintessence field, which interacts with 
the other such that the quintessence field has an oscillatory behaviour at the beginning of the 
matter dominated epoch and grows later so as to dominate the dynamics of the universe. None 
of the quintessence fields already there in literature has a proper physical motivation and it 
will be even more embarrassing to introduce a second field without any underlying physics. 
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So, the better arena is provided by a scalar tensor theory, such as Brans-Dicke theory, where 
one scalar field is already there in the purview of the theory and is not put in by hand. 
The relevant action in Brans-Dicke theory is given by 



S 



^d^x , (1.105) 



where \& is the Brans-Dicke scalar field, R is the Ricci scalar, Gq is the Newtonian constant 
of gravitation, uj is the dimensionless Brans-Dicke parameter, is the quintessence scalar 
field and f/(\l/,0) is the potential via which the two fields interact amongst themselves. We 
choose the form of 0) as V^(0)\1'^'^ which was used by Banerjee and Ram [52] for finding 
a solution of the graceful exit from infiationary paradigm in Brans-Dicke theory. With a 
slow roll approximation {<p^,(f) ~ 0), the conditions for an initially oscillating which grows 
during later stages are found out for two examples - a power law expansion and an exponential 
expansion of the scale factor a. 

For a power law expansion of the scale factor of the form 

a = aot", n > 1, (1.106) 

the condition on the potential V{(j)) comes out as 

iinVri^^-jt^-' (1.107) 

where A is a constant involving the parameters of the model. The subscript 'z' stands for 
some initial time {t = ti). 

Similarly, for the exponential expansion of the form 

a = aoe"*, ao,a>0, (1.108) 
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the condition on the potential comes out as 

(ln\/)"|.^0. (1.109) 

This indicates that for an exponentially expanding universe, the quintessence potential V{(f)) 
should be exponential at least in the beginning. 

In this paper two simple examples have been considered. However, this work can be ex- 
tended for more complicated kind of accelerated expansion of the universe. A key feature of 
this work is that here the numerical value of the Brans-Dicke parameter lj is not severely con- 
strained to very low values in order to drive the acceleration. For example, in the exponential 
expansion case, u is completely arbitrary whereas in power law expansion case it has been 
found that value of u can be determined in terms of other free parameters of the model. So, 
in both the cases u can be adjusted to some high value which is compatible with the limits 
imposed by solar system experiments |114j . 

In the second paper in this chapter entitled ^^An interacting scalar field and the recent 
cosmic acceleration" , the possibility of transition from deceleration to acceleration has been 
studied in the framework of a generalised scalar tensor theory. 

The dark matter and dark energy components are usually considered to be non-interacting 
and their evolutions are considered to be independent of each other. Recently Zimdahl and 
Pavon |115l 1116] showed that the interaction between the dark matter and dark energy com- 
ponents could be useful in solving the coincidence problem and an interacting scenario may 
provide a more general and better framework for obtaining an accelerating universe (see also 

[nzimiiin]) ■ 

In this paper, we introduce an interaction between the dark matter and the geometrical 
scalar field 0. The Brans-Dicke field equations are written in so called "Einstein frame" 
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as discussed earlier. The field equations in this version look tractable, although one has to 
sacrifice the equivalence principle as the rest mass of a test particle becomes a function of the 
scalar field (p. However, in this work, the metric has been transformed back to the original 
atomic units and the conclusions are drawn only from this. 
In place of the usual matter conservation equation 

p'm + 3Hpm = 0, (1.110) 

we have chosen an interaction between the dark matter and the scalar field (p of the form 

where pm is the matter density of the universe, H is the Hubble parameter, a is a positive 
constant. The overhead bars indicate that the equations are written in the conformally trans- 
formed version (Einstein's frame) where the transformed metric components are related to the 
original ones as 

9ti,y = <p9tMU ■ (1.112) 

The negative sign in equation fll.llip indicates that the energy is transferred from the dark 
matter component to the scalar field (p. It is found that this particular choice of interaction 
gives rise to a simple power law solution of the scale factor a - which can be either ever 
accelerating or ever decelerating depending upon the choice of parameters. Obviously we are 
not interested in such a scenario which does not provide the transition from deceleration to 
acceleration. One way out of this problem is to consider a generalization of Brans-Dicke theory 
where the parameter a; is a function of the scalar field (p rather than being a constant [57] . 
We make a choice of uj as 
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a being a positive constant. 

With this choice we obtain the solution for the scale factor in original units as 

AtT^ 

a = 7-. , 1-114 

(1 - 00^) ^ ' 

00, ^ are constants. 

From equation fll.ll4p one can easily find out the expression for q in the original version 



g = -i + ^2 ! ^ , 7.1. • (1-115) 



as 

[3^(l-0ot) + 0otP ■ 

The q vs. t plot shows the transition from positive value to a negative value for small values 
of a and also the nature of the curve is not critically sensitive to small changes in the value 
of a. 

However, from equation (11.114^ it is evident that at t — a, if and q all blow up. So 
the model is not valid upto infinite future and has a future singularity. So the model mimics 
a phantom model and has a big rip. 

Also, it is evident from equation f ll.llSp that {2uj + 3) has to be positive definite to allow a 
consistent model. But uj does not have any stringent limit and thus it is possible to adjust u 
to some high value, compatible with the limits imposed by the solar system experiments |114j . 
Also equations (11.1131) and fll.ll4p reveal that in the oo limit, a t^l"^ for small values 
of a. This is consistent with the notion that Brans-Dicke theory yields General Relativity in 
the infinite uj limit. Also for this model we have calculated the statefinder parameters {r, s} 
introduced by Sahni et. al [11], 1120] and found that the nature of the r vs. s plot for small 
values of a is similar to the one expected for scalar field quintessence models. 

The salient feature of this model is that no dark energy component is required here. The 
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interaction between the geometrical scalar field and matter component can drive the present 
acceleration of the universe. Although a particular form of interaction has been chosen, this 
is by no means unique and other forms of interaction can also be tried. 

The chapter 5 consists of one paper entitled Curvature driven acceleration : a utopia or 
a realityl" . This work explores the possibility whether geometry itself can serve the purpose 
of explaining the present accelerated expansion without having to resort to some exotic scalar 
field models. A substantial amount of work has been done along this direction which can 
provide an accelerating model for the universe [221 E31 El El El El ESI [ID] as discussed earlier. 
However, none of them could provide a smooth transition from deceleration to acceleration 
demanded by both theory }12H I122J and observations [123] . 

In the present work, we have written the field equations with a general f{R) and studied 
the model for two specific cases, namely, f{R) = i? — ^ and f{R) = e~^/^ . In both the cases 
it is shown that the required transition from deceleration to acceleration can be obtained. 
The relevant action for a general f{R) is 





For a spatially flat Robertson- Walker spacetime, where 



ds^ = dt^ - a^{t)[dr^ + r^dO^ + sin^ Odcp'^] , 




we write down the field equations for a vacuum universe taking L, 



m 



0. 



In the first example we make a choice of f{R) as 



f{R) = R 





where /i is a constant. This is exactly the form used by Carroll [5] and Vollick [6]. 
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For this particular choice, from the field equations we arrive at a relation 

1 



2H 



where the Ricci scalar R is given by 



1.119) 



(a a?\ , , 

i? = -6 - + — . 1.120 
\a j 

Equation fll.llQp is a highly nonlinear equation involving fourth order derivative of the scale 
factor a and it is difficult to obtain a complete analytic solution for a. But as we are inter- 
ested in the evolution of the deceleration parameter g, we translate equation (11.1191) into the 
evolution equation for q by using the relation 

a/a H 

Also, the time derivatives are replaced by the derivatives w.r.t. H . The resulting equation 
takes the form 

^gtt(^2_i)^2_2^2^t2(^+2)-igt^(g-l)(4g+7)-(2g+l)(g-l)2+if^(g-l)^ = , (1.122) 
where a f sign indicates differentiation w.r.t. H . 

This equation is also highly nonlinear and cannot be solved analytically. But if one can provide 
two initial conditions for q and q^ ^ a numerical solution is on cards. So, we pick up sets of 
values for q and q^ for H = 1 (i.e, the present values) from observationally consistent regions 
[TTl 1120] and plot q versus H numerically. The plot definitely shows a signature fiip in q from 
positive towards negative. Also, the nature of the curve is not critically sensitive to the initial 
conditions chosen. 

In the second example, the form of f{R) is chosen as 

f(R) = . (1.123) 
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Following the same method as before, the evolution equation for q is written as a function 
of H. With similar initial conditions for q and q^ at H = 1, q is numerically plotted against 
H and it is found that the curve has features similar to the previous example, i.e, q has a 
signature change from positive to negative in the recent past. The added feature of this latter 
example is that q has another signature flip from negative to positive direction indicating a 
decelerated expansion again in near future. 

So, in the present work both the examples indicate that gravity in its own right can lead to 
a late surge of accelerated expansion having a past deceleration which is essential for explaining 
nucleosynthesis and the structure formation of the universe. An added advantage of the second 
example is that in this case the universe re-enters a decelerated expansion phase in the near 
future and thus ^phantom menace^ is avoided, i.e, the universe does not attain infinite rate of 
expansion in a finite future. 

There is already some criticism of gravity regarding its stability [TT], but still there 
are reasons to be optimistic about a curvature driven acceleration. It deserves mention that 
although the term has been quite widely used in the literature, it suffers from the drawback 
that at i? = 0, i.e, in the late stage of the evolution, the model has a singularity. The other 
example, i.e, e~^^^ is regular everywhere. The present work deals with a vacuum universe, 
but in a more general framework one has to either put in matter or should derive the relevant 
matter from curvature itself. Thus the modified theory of gravity provides a platform for more 
ambitious and detailed work in future. 
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Chapter 2 

A Simple Quintessence Field 
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2.1 Acceleration of the Universe with a Simple Trigono- 
metric Potential 

(Journal reference : N. Banerjee and S. Das, Gen. 
Rel. Grav., 37, 1695 (2005); astro-ph/0505121, ) 
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2.1.1 Introduction 

Over the last few years, there are growing evidences in favour of the scenario that the universe 
at present is expanding with an acceleration. The supernovae project \T] and also the Max- 
ima [2j and Boomerang [3] data on cosmic microwave background (CMB) strongly suggest 
this acceleration. The very recent WMAP data [1] also seem to confirm this. This result 
is indeed counter-intuitive, as gravity holds matter together and it might be expected that 
in the absence of any exotic field, the universe should be decelerating. As this acceleration 
could be brought about by an effective negative pressure, the first choice of candidate for this 
'dark energy' had been the 'cosmological constant' or a time varying cosmological parameter 
A(t). Due to well-known reasons, A has fallen from grace (for an excellent uptodate review, 
see [5],[6])- A scalar field with a positive definite potential can indeed give rise to an effec- 
tive negative pressure if the potential term dominates over the kinetic term. The pressure 
to density ratio for the scalar field, as required by the supernovae observations, is given as 
< — I (See ref [7] and references therein). This source of energy is called the quintessence 
matter (Q-matter). In this context, non-minimally coupled scalar fields had been investigated 
thoroughly to check if they could drive an accelerated expansion [8]. Brans-Dicke's scalar field 
appears to generate sufficient acceleration in the matter era, but it has its problems in the 
earlier evolution [9]. A viscous fluid along with a Q-matter could also be a useful candidate 
and this appears to solve the coincidence problem also [10] . This coincidence problem, i.e, 
why the Q-matter dominates only recently, was tackled in the so called tracker solutions pT] 
where the scalar field energy density runs parallel to the matter energy density from below 
through the evolution and gets to dominate only during later stages. Very recently Chaplygin 
gas, which has a nonlinear contribution of the energy density to the dynamics of the model. 
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has also been invoked [12]. Most of these models do exhibit an accelerated expansion in the 
matter-dominated regime. 



It deserves mention that the same model should have a deceleration in the early phase of 
matter era in order to provide a perfect ambience for structure formation. Furthermore, the 
accelerated phase is perhaps only a very recent one. There are observational evidences too 
that beyond a certain value of the redshift {z ~ 1.7), our universe had been going through a 
decelerated expansion [13]. This indication is indeed reassuring, as the formation of structure 
in the universe is better supported by a decelerating model. This is because local inhomo- 
geneities will grow and become stable from the seeds of density fluctuation only if the force 
field is attractive. 

Amendola [H] has argued that all the required structure formation and other relevant 
observations regarding the supernovae could well be explained even if the alleged acceleration 
of universe started quite a long time back, even beyond z = 5. However, this work also shows 
that the model requires both an accelerated and a decelerated phase of expansion. But a more 
recent work by Padmanabhan and Roychowdhury [15] shows a striking result. It indicates that 
if we take the complete data set, i.e, acceleration upto a certain z and deceleration beyond 
that (i.e, for higher z), then only this conclusion of the change of signature of the deceleration 
parameter holds. On the other hand, the individual data sets of the high and low redshift 
supernovae may well be consistent with a decelerating universe without any 'dark energy'. 

So indeed we are in need of some form of fields which governs the dynamics in such a 
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way that the deceleration parameter becomes negative well into the matter era. One such 
Q-matter had been given by Sen and Sethi [16] where they include a potential which is a 
'double exponential' of the scalar field. They obtained a scale factor which is a sine hyperbolic 
function of time in the matter dominated regime. The deceleration parameter (g) indeed has 
a sign fiip and with a little fine-tuning, the scale factor can grow during the early stages as 
t^/^ which is indeed the usual solution for the Einstein equations for a fiat FRW spacetime for 
pressureless dust. 

In the present work, we adopt the following strategy. We choose a form of g as a function of 
the scale factor a so that it has the desired property of a signature fiip. Then with this input, 
the scalar field and the required potential are found out. It turns out that a fairly simple 
trigonometric potential does the needful. The origin of the scalar potential, however, cannot 
be indicated. Surely this is not the ideal way to find out the dynamics of the universe, as here 
the dynamics is assumed and then the fields are found out without any reference to the origin 
of the field. But in the absence of more rigorous ways, this kind of investigations collectively 
might finally indicate towards the path where one really has to search. This 'reverse' way of 
investigations had earlier been used extensively by Ellis and Madsen [17] for finding out the 
potential driving infiation, i.e, an accelerated phase of the universe at a very early stage of its 
evolution. 



2.1.2 Results 

For a spatially fiat Robertson- Walker spacetime 



(2.1) 
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the deceleration parameter q is given by 

(2.2) 

where a is the scale factor of the universe and is a function of the cosmic time 't' alone. 
In order to get a model consistent with observations, one needs an expanding universe, i.e, a 
positive Hubble parameter H = ^ throughout the evolution, but a deceleration parameter q, 
unlike being a positive constant throughout the matter era at g = 1/2 as believed until the 
recent observations, should be a function of the scale factor ( or that of t ). Furthermore, this 
functional dependence should be such that q undergoes a transition from its positive phase to 
a negative one in the matter dominated period itself. It is thus imperative that the scale factor 
cannot have a simple power- law behaviour. If a ~ t", the universe will have an accelerated or 
a decelerated expansion for n > 1 or n < 1 respectively throughout the period. 

In the quest for a varying q consistent with observations, in the same line as that floated 
by Ellis and Madsen, we propose the relation 

where p is a constant. It is found that for a certain range of negative values of p, this works 
remarkably well. 

The equation (12.31) integrates to yield 

H=- = A(1 + an (2.4) 
a 

where A is an arbitrary constant of integration. A is taken to be positive, which ensures the 
positivity of the Hubble parameter (the expansion of the universe is never denied!) irrespective 
of the signature or value of the constant p. 
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It is found that for values of p between -2 and -1, the model shows exactly the behaviour 
which is desired (as shown in Figure EH]) . 




■ " ■ >H 



Figure 2.1: Plot of q vs. z for spatially fiat dust dominated R-W model. 

In what follows, we work out the problem completely for p = —3/2, for which the model 
works with a non minimally coupled scalar field with a potential expressed as a simple trigono- 
metric function of the scalar field. 

As the interest is in a matter dominated universe, the fluid is taken in the form of a 
pressureless dust. The Einstein equations for the space-time given by equation (12. ip are, 

3^ = P + ^0' + m, (2.5) 

2^ + ^ = -^0^ + n0) (2.6) 
a 2 

where p is the density of matter, is the scalar field and V^(0) is a scalar potential. 
The wave equation for the scalar field is 

+ 3-0 + K'(</)) =0 (2.7) 
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In all equations an overhead dot implies differentiation w.r.t. time and a prime is that w.r.t. 
the scalar field (p. The matter conservation equation, which can in fact be obtained from these 
three equations in view of the Bianchi identity, yields 

P = ^, (2.8) 

Po being a constant. So we have three equations to solve for four unknowns. 

We assume the deceleration parameter as given in equation (12.31) . which can be integrated 
twice to give if = d/a as in equation (12.41) and the scale factor as 

a = [e-^P* - (2.9) 



Now, the system of equations (12.51) . (12.61) and (12.71) is closed with the assumption of equation 
(12. 3p or equivalently equation (12. 4p . So in order to solve the system completely, the parameter 
'p' should have a fixed value. We choose p = — |, as it yields q = 0.5 for a very low value of 
where oq is the present value of the scale factor. The physical motivation for choosing this 
value of q for early matter dominated epoch is that for a spatially fiat FRW model with p = 
without any Q-matter indeed has g = 0.5 and that the transition from radiation to matter 
dominated epoch for this value of q is well studied p8]. 

With p = — |, equations (12.51) and (12. 6p are used to eliminate V{(j)), and can be 
calculated to be 

^ = V3Aa^ = -^^ (2.10) 
[e 2 — IJ 2 

The scalar field is found out by integrating equation (I2.10p as, 

0= ^tan-i(e+3^*/2-l)5, (2.11) 
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The potential V{(f)) can also be calculated from equations fl2.5p and (12. 6p first as a function 
of time and by the use of equation (12. lip as a function of (p as, 

V{^) = ^-^cot\y^) + 3A^ (2.12) 

Now one has the complete set of the solutions, a = a{t), = (p{t), p = p{t) and V = V{(f)) 
for p = —3/2. The solutions, when plugged in the field equations, namely (12.50 . (12. 6p and 
(12. 7p . satisfy all of them provided 

po = 3^2. (2.13) 

From equations (12. 5p and (12. 6p . we note that the contribution from the quintessence field 
(f) towards the density and effective pressure are given as 

P0 = ^0' + n0), (2.14) 

and 

P<p = - V{<P) (2.15) 

respectively. From these, one can write down the expressions for the dimensionless density 
parameters flm = ff^ and = respectively for the visible matter and the Q-matter. 

With p = —3/2, the present model yields 

[l + (l + z)-3/2]2 ' ^^-^^^ 

and 



1 ) 



(2.17) 
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where z is the redshift parameter given by 
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1 I "0 
1 + z = — , 

a 



(2.18) 



ao being the present value of the scale factor, fi^o, the present value of Qm, comes out to be 
0.25 and fl^g = 0.75. These values are well within the constraints of 0.2 < flmo < 0.8 [T9].[5]. 
Figure 12.21 shows that increases with z, i.e, decreases with the evolution of the universe. 

starts dominating over Qm roughly a.t z = 0.8. At the earlier epoch, i.e, at high z, is 
very small, allowing a conducive matching onto the radiation era for the perfect ambience for 
nucleosynthesis. 
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Figure 2.2: Plot of vs. z. 



Unlike Newtonian gravity, general relativity ensures that the pressure also contributes in 
driving the acceleration of the model. From equations (12.51) and (12. 6p . one has 



d 1 
a D 



(2.19) 
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where pt and pt are the total effective density and pressure respectively. If the pressure is 
connected with the density by the relation 



p = wp, 



(2.20) 



the model will accelerate ( a > ) only if 



Pm+Pcj, 1 

Wt = < — . 

Pm+ Ptf, 3 



(2.21) 



The subscripts 't', 'm' and '0' stand for total, normal fluid distribution and the Q-matter 
respectively. For a matter dominated universe, Pm = and hence Wm = 0. This particular 



model gives 



+ Pr, 



1 + (1 + ^)^ 

-1 -2(l + z)i 

-[l + (l + z)t]-l 



(2.22) 
(2.23) 



The evolution of Wt versus z and w^p versus z are shown in figure 12. 3[ which indicates that Wt 
attains the required value of — | or less only close to z = 0.5. Beyond that, Wt is less negative, 
and the universe still decelerates. 




O.i 



1.5 




Figure 2.3: Plot of (a) Wt vs. z and (h) w^i) vs. z. 



The value of w^q, i.e, the value of the equation of state parameter for the scalar field at 
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2; = as given by the present model and as indicated by figure 12.3( b) is definitely within the 
constraint range [T9] . 
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Figure 2.4: Plot of vs. z for spatially flat dust dominated R-W model. 

We also plot the rate of change of Wt against z ( as shown in figure 12.41 ) . It shows that 
^ is still negative at the present epoch, but the magnitude of ^ is decreasing. 

The solution for the scale factor is good enough to allow the density contrast to grow 
favourably for the formation of large scale structure. Figure [231 shows the growth of linearized 
density perturbation in this model and evidently indicates that it grows linearly with the scale 
factor during later stages as expected for the matter dominated epoch [18]. 

In the absence of the final form of the quintessence matter, search for the relevant form of 
potential will continue and the present investigation is one of them. In view of the high degree 
of non-linearity of Einstein's equations, exact solutions always play a vital role as piecewise 
solutions have the problem of proper matching at different interfaces. The present model 
shows that inspite of the severe constraints imposed by observations, one can still find an 
exact FRW model which gives values for the relevant parameters like g, w, VL^, Vtm etc. safely 
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Figure 2.5: Plot of density contrast D vs. a where D = 

within the range given by observations. It also has the merit of having a single analytical 
expression for q = q{a), which gracefully transits from its positive phase to the negative one 
and adds to the list of quintessence potentials that serve the purpose of modelling a presently 
accelerating universe [20]. Definitely the model has problems, particularly that of fine tuning, 
but infact, all quintessence models have some such problems. 
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3.1 Spintessence : A Possible Candidate as a Driver of 
the Late Time Cosmic Acceleration 
(Journal reference : N. Banerjee and S. Das, Astro- 



phys. Space Sci., 305, 25 (2006); gr-qc/0512036,) 
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3.1.1 Introduction 

The stunning results of the observations on the luminosity - redshift relation of some distant 
supernovae \T^,^, that the universe is currently undergoing an accelerated phase of expansion, 
poses a serious challenge for the standard big bang cosmology. As the standard gravitating 
matter gives rise to an attractive field only, this challenge is negotiated in the standard model 
by invoking some field which gives rise to an effective negative pressure. For many a reason, 
a dynamical 'dark energy' is favoured against the apparently obvious choice of a cosmological 
constant A for providing this negative pressure [3]. This dynamical dark energy is called a 
quintessence matter. It cannot be overlooked that a successful explanation of the formation of 
structures in the early matter dominated era crucially requires an effectively attractive gravi- 
tational field and thus a decelerated phase of expansion must have been witnessed by an early 
epoch of matter dominated universe itself. Very recently Padmanabhan and Roy Chowdhury 
[1] showed that the data set, only showing the accelerated phase of expansion, can well be 
interpreted in terms of a decelerated expansion in disguise. The acceleration only becomes 
meaningful if the full data set shows deceleration upto a certain age of the universe and an 
acceleration after that (see also ref [5j). In keeping with such theoretical requirements, actual 
observations indeed indicate such a shift in the mode of expansion - deceleration upto a higher 
redshift regime (about z ~ 1.5) and acceleration in more recent era, i.e, for lesser values of z [H]. 

This observation indeed came as a relief, as the formation of galaxies could proceed un- 
hindered in the decelerated expansion phase. It also requires all quintessence models to pass 
through certain fitness tests, such as the model should exhibit a signature flip of the decel- 
eration parameter q from a positive to a negative value in the matter dominated era itself. 
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Quite a few quintessence models do exhibit such a signature flip of the deceleration parame- 
ter. One very attractive model was that of 'spintessence' proposed by Boyle, Caldwell and 
Kamionkowski [7]. It essentially works with a complex scalar field and a potential, which is a 
function of the norm of the scalar field. The scalar field, 

%jj = (p^ + i(p2, (3.1) 

can be written as 



= 06*"*, (3.2) 

i.e, the complex part is taken care of by a phase term. This kind of a scalar field is already 
known in the literature in describing a 'cosmic string' [8] although it has to be noted that a 
cosmic string has a very specific form of the potential V{(j)), whereas, the relevant form has to 
be found out for a quintessence model. Boyle et al discussed the so called spintessence model 
in two limits separately, namely for a high redshift region and also for a very low redshift 
region. Apparently this model gives exactly what is required, an acceleration for the low z 
limit whereas a deceleration for a high redshift limit. In the present investigation, it is shown 
completely analytically that the model indeed works. For high value of the scale factor a, 
the deceleration parameter q is negative whereas for a low value of a, q is positive. The most 
important feature is that the value of a, where the sign flip of q takes place, can be analytically 
expressed in terms of the parameters of the theory and constants of integration. 
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3.1.2 Field Equations and Results 

If we take the scalar field as given by equation (13.21) . Einstein's field equations become 

3^=P + U' + Ico'^' + Vi<P) , (3.3) 

2^ + ^ = -l0^-l^V + m, (3.4) 
a 2 2 

where a is the scale factor, p is the energy density of matter, V is the scalar potential which 
is a function of amplitude of the scalar field, and the phase u is taken to be a constant. An 
overhead dot implies differentiation w.r.t. the cosmic time t. In a more general case, u might 
have been a function of time. 

The matter distribution is taken in the form of dust where the thermodynamic pressure p 
is equal to zero. This is consistent with the 'matter dominated' epoch. This leads to the first 
integral of the matter conservation equation as 

P = ^. (3.5) 
Variation of the relevant action with respect to (f) yields the wave equation 

^ + 3-<j) + V'{(f)) = u^(j), (3.6) 
a 

where a prime is a differentiation w.r.t. (p. 
As the scalar field actually has two components, a third conservation equation is found as 

Aq being a constant. 

In the spintessence model described by Boyle et al [7] , it has been considered that u; is a very 
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slowly varying function of time which indeed can be considered as a constant over the entire 
period of dust dominated era. So for constant uj, equation ( 13 .7^ yields 

0^0^ = A , (3.8) 

A being a positive constant. 

It deserves mention that only two equations amongst fl3.5p . fl3.6p and fl3.8p are independent 
as any one of them can be derived from the Einstein field equations with the help of the other 
two in view of the Bianchi identities. So, we have four unknowns, namely, a, p, and V{(f)) 
and four equations, e.g, (13. 3p . fl3.4p and two from (I3.5l) - fl3.8p . So this is a determined problem 
and an exact solution is on cards without any input. From equations (13.31) . (13.41) and (13. 5p 
one can write 

(3.9) 



a po 1 -2 1 



which takes the form 



a^'^ = -l-niH\ (3.10) 

where H = ^, the Hubble parameter, 

I = |(Po + io'^A) and m = ^ are positive constants. 

In deriving this, the equation (13.80 has been used to eliminate (f) and 4> in terms of a and 
a. Now we make a transformation of time coordinate by the equation 

As ^ is positive definite ( a is the scale factor and cannot take negative values ), we find 
that a; is a monotonically increasing function of t. So one can use x as the new cosmic time 
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without any loss of generality and deformation of the description of events. 
In terms of x, equation ( I3.10p can be written as 

= -l- mH^ , 3.12 

ax 

which can be readily integrated to yield 

H = n tan{(3 — mnx) , (3.13) 

where = ^, is a positive constant, and /5 is a constant of integration. 

Now we use this if as a function of the new cosmic time variable x and find out the 
behaviour of the deceleration parameter q, which is defined as 

H W 

where a dagger indicates differentiation w.r.t. x. 
Equation fl3.12p now yields 

g = -l + ^(^ + m) . (3.15) 
Both m and / are positive constants. So q has a 'zero', when 

ai^ = ^ + m. (3.16) 

The suffix 1 indicates the values of the quantities for g = 0. Furthermore, equation (13.151) can 
be differentiated to yield 

rln if/ 1 

(3.17) 



dq 
da 



1 

af 



^ + 3m 



at the point q = 0, a = ai and H = Hi. In deriving the equation (13.171) . equations (I3.10p 
and fl3.16p have been used. The last equation clearly shows that g is a decreasing function of 
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a atleast when g = 0. So g definitely enters a negative value regime from a positive value at 
a = ai and H = Hi. 

For the sake of completeness, the solution for the scale factor a can be found out by 
integrating equation (13.131) as 



m 



mnx) I 



1/3 

(3.18) 



3.1.3 Discussion 



So evidently, the spintessence model proposed by Boyle et al [7J passes the 'fitness test', the 
deceleration parameter enters into a negative value in a "finite past" . In view of the high degree 
of non linearity of Einstein equations, exact analytic solutions are indeed more dependable, 
and the present investigation provides that in support of a spintessence model. Also, the 
constants of the theory ( such as uj ) and the constants of integration ( such as / and m ) are 
still free parameters and hence provides the 'comfort zone' for fitting into the observational 
results. Another feature of this study is that the results obtained are completely independent 
of the choice of potential V = V {(/)). This feature provides a bonus, as different potentials, used 
as the quintessence matter, are hardly well-motivated and do not have any proper physical 
background. Recently quite a few investigations show that a complex scalar field indeed 
serves the purpose of driving a late time acceleration [9J. But most of these investigations 
either invoke the solution in some limit ( such as for a large a ) or use some tuning of the 
form of the potential. The present investigation provides a better footing for them as it shows 
the transition of q analytically. It also deserves mention that Bento, Bertolami and Sen [10] 
showed the effectiveness of a Chaplygin gas as a quintessence matter. It is interesting to note 
that under some assumptions this kind of a fluid formally resembles a complex scalar field as 
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discussed in this work. 
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4.1 A Late Time Acceleration of the Universe with Two 
Scalar Fields : Many Possibilities 
(Journal reference : N. Banerjee and S. Das, Mod. 
Phys. Lett. A, 21, 2663 (2006); gr-qc/0605110, ) 
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4.1.1 Introduction 

The present cosmic acceleration is now generally believed to be a certainty rather than a 
speculation. The recent data on supernovae of type la suggested this possibility quite strongly 
[H [21 El H] and the most trusted cosmological observations, namely that on Cosmic Microwave 
Background Radiation [5], El [71 El [9], appear to be quite compatible with an accelerated 
expansion of the present universe. The natural outcome of these observations is indeed a 
vigorous search for the form of matter which can give rise to such an expansion, as a normal 
matter distribution gives rise to an attractive gravity leading to a decelerated expansion. This 
particular form of matter, now popularly referred to as "dark energy" , is shown to account 
for as much as 70% of the present energy of the universe. This is also confirmed by the highly 
accurate Wilkinson Microwave Anisotropy Probe (WMAP) [TOl [Til [121 [131 [H] . A large number 
of possible candidates suitable as a dark energy component have appeared in the literature. 
Excellent reviews on this topic are available [TBI [T7] . Most of the dark energy candidates 
are constructed so as to generate an effective pressure which is sufficiently negative driving an 
accelerated expansion. The alleged acceleration can only be a very recent phenomenon and 
must have set in during the late stages of the matter dominated expansion of the universe. 
This requirement is crucial for the successful nucleosynthesis in the radiation dominated era 
as well as for a perfect ambience for the formation of structure during the matter dominated 
era. Fortunately, the observational evidences are also strongly in favour of a scenario in which 
the expansion of the universe had been decelerated (the deceleration parameter g > 0) for 
high redshifts and becomes accelerated {q < 0) for low values of the redshift z [18j. So the 
dark energy sector should have evolved in such a way that the consequent negative pressure 
has begun dictating terms only during a recent past. 
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This so easily reminds one about the inflationary universe models where an early accel- 
erated expansion was invoked so as to wash away the horizon, fine tuning and some other 
associated problems of the standard big bang cosmology. The legendary problem, that the 
inflationary models themselves had, was that of a "graceful exit" - how the accelerated expan- 
sion gives way to the more sedate (g < 0) expansion so that the universe could look like as we 
see it now. For a very comprehensive review, we refer to Coles and Lucchin |19| or Kolb and 
Turner [20]. The "graceful exit" problem actually stems from the fact that for the potentials 
driving inflation, the phase transition to the true vacuum is never complete in a sizeable part 
of the actual volume of the universe. The attempts to get out of this problem involved the in- 
troduction of a scalar fleld which slowly rolls down its potential so that there is sufficient time 
available for the transition of phase throughout the actual volume of the universe. The current 
accelerated expansion thus poses the problem somewhat complementary to the graceful exit - 
that of a "graceful entry" . 

The present work addresses this problem, and perhaps provides some clue regarding the 
solution of the problem. The basic motivation stems from existing literature on infiation. 
Mazenko, Wald and Unruh [2T] showed that a classical slow roll is in fact invalid when the 
single scalar field driving infiation is self interacting. It was also shown that a slow roll with a 
single scalar field puts generic restrictions on the potentials driving infiation |22j . This kind of 
problems led to the belief that for a successful infiationary model, one needs to have two scalar 
fields [23]. The present work uses this idea of utilising two scalar fields, one of them being 
responsible for the present acceleration of the universe and is called the quintessence field. 
The second one interacts nonminimally with the former so that the quintessence field has an 
oscillatory behaviour at the early matter dominated epoch but indeed grows later to dominate 
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the dynamics of the more recent stages of the evolution. If such a behaviour is achieved, some 
clue towards the resolution of the graceful entry problem or the coincidence problem may be 
obtained. There are quite a few quintessence potentials already in the literature [I5] which 
drives a late time acceleration, but none of them really has an underlying physics explaining 
their genesis. As one is already hard pressed to find a proper physical background of the 
quintessence field, it will be even more embarassing to choose a second field without any 
physical motivation. Naturally, the best arena is provided by a scalar tensor theory, such as 
Brans-Dicke theory, where one scalar field is already there in the purview of the theory and is 
not put in by hand. It deserves mention that the Brans-Dicke scalar field was effectively used 
in "extended infiation" in order to get a sufficient slow roll of the scalar field f2M [25] . Later 
Brans-Dicke theory was used for finding a solution of the graceful exit problem with a large 
number of potentials [26] , where the infiaton field evolves to an oscillatory phase during later 
stages. 

In the next section we write down the field equations in Brans-Dicke (B-D) theory with 
a quintessence field 0, the potential V{(j)) driving acceleration being modulated by the B-D 
scalar field as With a slow roll approximation, the conditions for an initially 

oscillating 0, which grows only during later stages, are found out for two examples, a power 
law expansion and an exponential expansion of the scale factor. The particular form of V{(f)) 
is quite irrelevant in this context, the conditions only put some restrictions on the constants 
of the theory and the parameters of the model. So the form of the potential is arbitrary to 
start with, only the conditions on the parameters and the 'value' of V{(f)) has to be satisfied 
and hence many possibilities are opened up to accommodate a physically viable potential as 
the driver of the late time acceleration. However, in some cases, this could restrict the form 
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of V{(f)) as well. In the last section, we make some remarks on the results obtained. 
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4.1.2 A Model with a Graceful Entry 

The relevant action in Brans-Dicke theory is given by 

S^J [j^^ - C.^^ - - Uiij, 0) + Lrn]V^d^X (4.1) 

where Gq is the Newtonian constant of gravitation, uj is the dimensionless Brans-Dicke param- 
eter, R is the Ricci scalar, ■0 and cf) are the Brans-Dicke scalar field and quintessence scalar 
fields respectively. If we now choose U{^, (p) as V{(l))^~'^ as explained, the field equations, in 
units where SuGq = 1 , can be written as 

3H^ + V'{(l))ijj-^ = , (4.3) 

{2uj + 3) (V^ + 3HiP) = (/3 - A)V{(l))ip-^ + p , (4.4) 

where a dot represents differentiation with respect to time t and a prime represents differenti- 
ation with respect to the scalar field 4>. As we require the potential U (p) to grow with time 
so that the effective negative pressure dominates at a later stage, P should be negative for a i/j 
growing with time or positive for a ip decaying with time. The field equations are written in 
the slow roll approximation, i.e, where 0^ and are neglected in comparison to others. H — ^ 
is obviously the Hubble parameter. As we dropped the field equation containing stresses, we 
can use the matter conservation equation as the fourth independent equation which yields on 
integration 
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Po being a constant. This is so as the fluid pressure is taken to be zero as we are interested in 
the matter dominated era. 

Using the expression for V{(f))ip~^ from equation (14. 2 p in equation (14.41) . we can write 

(2. + 3)|| + 3^] ^ - + SHt - I J, ^ (, - 5)^ . (4.6) 

If the scale factor a is known, this equation can be integrated to yield the Brans-Dicke field 

Equations fOj) and fl4i|) yield 



(4.7) 



Hence, if we define 



then 



/(0o)-/(0.)= r F{t)dt, (4.9) 
Jti 

where 

_ _(2. + 3)(^ + 3g^)-p 

3(/5-4)/7 ^ ' 

and subscripts 'o' and 'i' stands for the present value and some initial value, such as the onset 
of the matter dominated phase of evolution. 

For a given a = a{t), therefore, equation (14. 6 p can be used to find ip, which in turn, with 
equations (14.81) and (14.91) determines /(0). Now, these equations can be used to put bounds 
on the values of derivatives of the potential, which would ensure that the dark energy has an 
oscillating phase in the early stages. 

The complete wave equation for the quintessence field is 
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+ 3i/0 + r(0)V^-^ = 0. 

The condition for a small oscillation of (f) about a mean value is V'{(j)) = 0. This provides a 
kind of plateau for the potential which hardly grows with evolution and hence the dynamics 
of the universe is practically governed by the B-D field ip and the matter density p. We find 
the condition for such an oscillation of at some initial epoch by choosing V'{(f)) ~ which 
yields 




During later stages, the universe evolves according to the equations (14.21) - (14.41) where V'{(j)) ^ 
0, and the scalar field slowly rolls along the potential so that the quintessence field takes an 
active role in the dynamics and gives an accelerated expansion of the universe. 

Two examples, one for power law and the other exponential expansion, will be discussed 
in the present work. 



I. Power law expansion :- 

If a = aot" where n > 1, the universe expands with a steady acceleration, i.e, with a constant 
negative deceleration parameter q = — ■ 
With this. 



H=^- = '}, (4.11) 

a t 



and equation (14. 6 p has the form 



Cj's being constants given by, 
ci = (2cj+3), C2 = 3n{2uj-(3+7), cg = ^(/5-4), = -3n\f3-A), cs = (/3-5)^ . (4.13) 
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The simplest solution for ip in equation (14.121) is 

^ = ^ot'-'" , (4.14) 

ipo being a constant. 

The consistency condition for this is, 

(2cu+3)(2-3r2)(l-3n)+3n(2tu-/3+7)(2-3n) + ^(/3-4)(2-3r2)2-3n2(/5-4) + (/5-5) ^° - ^ 



From equations (14.91) and (I4.10p . 

f{^o)-f{<P^) = Dit^~'--tl-'-), (4.15) 
where D is a constant involving q's, i.e, n, u, ipo, Po etc. given by 



(2a; + 3)^o(2 - 3n) - 



3n(/3-4)(2-3n) 

From the form of potential V = V{(j)), f{(j)) can be found from the relation (14. 8p . So, the 
different constants will be related by equation (14.151) . 

If the quintessence field oscillates with small amplitude about the equilibrium at the be- 

i 1, which 

puts more constraints amongst different parameters of the theory. 
Equations (griD and fHTTD yield 



ginning, i.e, at t = and grows during the later stages of evolution, then 



= A {In vy , 



where A is a constant given by 



(2^ + 3)^o(2 - 3n) - ^ 



PO 



A = — ^ ^ = -D(2 - 3n) 

3n{p - 4) ^ ' 
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So, 

3H^ 3n 3n 

The condition for the small oscillation of 6 close to t = t,- is 



3H(f) 



which gives the condition on {In V)" as 



(4.16) 



II. Exponential expansion :- 

Similarly for an exponential expansion at the present epoch, the constraints can be derived. 
For such an expansion, 

a = aQ e , 

where Oq, a are all positive constants. 
Then, 

H=^ = a and p = ^e'^"* . 
a an 



Equation (14.61) with this solution has the form 



where 6j's are constants given by 



6i = (2cj+3), b2 = 3a{2iu-l3+7), 63 = ^(/3-4), h = -3a\(3-A), h = {f3-5)^ . {A.U 
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A simple solution for ip is 

= ^oe-'"* , (4.19) 
ipo being a constant. The consistency condition for this is, 

9a\(3 + - 2cu - 4) - 3a\(3 - 4) + (/3 - 5)-^ = . (4.20) 
2 a^ijjo 

Equations (14. 9 p and (14.101) will put restrictions on the parameters of potential by 

(4,21) 

The condition ^ —1 for the small oscillation of the scalar field at t = tj, with the help of 
equations (14.71) and (I4.19p . leads to the interesting result 

{In V)"\ . ^ , (4.22) 

which indicates that for an exponentially expanding present stage of evolution, the quintessence 
potential behaves exponentially at least at the beginning. 

Thus, if = V{(j)) is given, then equations (14.151) and (14.211) will help finding out the 
bounds on the values of the scalar field and the constants appearing in V{(f)) in terms of the 
initial and final epochs. 

The conditions (I4.15P or (14.211) put bounds on the potential so that the quintessence field 
has an oscillatory behaviour in the beginning. As the solutions in the examples considered 
has accelerated expansion, the Q-field has a steady growth at later stages. For a power law 
expansion, the growth has some arbitrariness as V{(j)) is not specified. For the exponential 
expansion, however, the growth of is governed by an exponential potential. 
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4.1.3 Conclusion 

It deserves mention that non-minimally coupled scalar fields were utilised by Salopek et al|27j 
and Spokoiny [28| [29] in the context of an early inflation, where the field had an oscillation 
or a plateau at the end of the inflation. In the present case, we need just the reverse for the 
Q-field, and that is aided by the non-minimally coupled field i/). 

We see that for a wide range of choice of V{(f)), a power law acceleration is on cards, only 
the value of V{(j)) at some initial stage is restricted by equation fl4.16p . For an exponential 
expansion of the scale factor, however, the potential V{(j)) has to be an exponential function 
of ( in view of equation (I4.22p ). This investigation can be extended for more complicated 
kinds of accelerated expansion. 

It is true that general relativity is by far the best theory of gravity and the present cal- 
culations are worked out in Brans-Dicke (B-D) theory, but this should give some idea about 
how a second scalar field may be conveniently used to get some desired results. Although B-D 
theory lost a part of its appeal as the most natural generalization of general relativity (GR) 
as the merger of B-D theory with GR for large u limit is shown to be somewhat restricted 
|30] . it still provides useful limits to the solution of the cosmological problems [24j. Another 
feature of the present work is that the numerical value of u required is not much restricted. 
For power law expansion, u is restricted by equation fl4.16l) which clearly shows that it can be 
adjusted by properly choosing values of some other quantities, whereas for exponential expan- 
sion, equation (14.221) shows that uj is arbitrary. This is encouraging as it might be possible to 
get an acceleration even with a high value of uj, compatible with local astronomical observa- 
tions ^31j. B-D theory had been shown to generate acceleration by itself [32], although it had 
problems with early universe dynamics. B-D theory with quintessence or some modifications 
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of the theory [331 EH ESI EH EZ] were shown to explain the present cosmic acceleration, but 
all these models, unlike the present work, required a very low value of uj, contrary to the local 
observations. 
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4.2 An Interacting Scalar Field and the Recent Cosmic 
Acceleration 

(Journal reference : S. Das and N. Banerjee, Gen. 
Rel. Grav., 38, 785 (2006); gr-qc/0507115, ) 
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4.2.1 Introduction 

Over the last few years, the speculation that our universe is undergoing an accelerated expan- 
sion has turned into a conviction. The recent observations regarding the luminosity - redshift 
relation of type la supernovae [H [2l [3|, H] and also the observations on Cosmic Microwave 
Background Radiation (CMBR) [SI El [71 El E] very strongly indicate this acceleration. These 
observations naturally lead to the search for some kind of matter field which would generate 
sufficient negative pressure to drive the present acceleration. Furthermore, observations re- 
veal that this unknown form of matter, popularly referred to as the "dark energy" , accounts 
for almost 70% of the present energy of the universe. This is confirmed by the very recent 
Wilkinson Microwave Anisotropy Probe (WMAP) data [101 [IH [121 [iSl [Tl| . A large number 
of possible candidates for this "dark energy" component has already been proposed and their 
behaviour have been studied extensively. There are excellent reviews on this topic [T6| [38]. 

It deserves mention that this alleged acceleration should only be a very recent phenomenon 
and the universe must have undergone a deceleration (deceleration parameter q = —-^j^ > 0) 
in the early phase of matter dominated era. This is crucial for the successful nucleosynthesis 
as well as for the structure formation of the universe. There are observational evidences too 
that beyond a certain value of the redshift z { z ^ 1.5 ), the universe surely had a decelerated 
phase of expansion [T8]. So, the dark energy component should have evolved in such a way 
that its effect on the dynamics of the universe is dominant only during later stages of the 
matter dominated epoch. A recent work by Padmanabhan and Roy Choudhury [391 [IQ] shows 
that in view of the error bars in the observations, this signature fiip in q is essential for the 
conclusion that the present universe is accelerating. 

So, we are very much in need of some form of a field as the candidate for dark energy, which 
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should govern the dynamics of the universe in such a way that the deceleration parameter q 
was positive in the early phases of the matter dominated era and becomes negative during 
the later stages of evolution. One of the favoured choices for the "dark energy" component 
is a scalar field called a quintessence field ( Q-field ) which slowly rolls down its potential 
such that the potential term dominates over the kinetic term and thus generates sufficient 
negative pressure for driving the acceleration. A large number of quintessence potentials 
have appeared in the literature and their behaviour have been studied extensively ( for a 
comprehensive review, see jTH] ). However, most of the quintessence potentials do not have 
a proper physical background explaining their genesis. In the absence of a proper theoretical 
plea for introducing a particular Q-field, non-minimally coupled scalar field theories become 
attractive for carrying out the possible role of the driver of the late time acceleration. The 
reason is simple; the required scalar field is already there in the purview of the theory and 
does not need to be put in by hand. Brans - Dicke theory is arguably the most natural 
choice as the scalar - tensor generalization of general relativity (GR) because of its simplicity 
and a possible reduction to GR in some limit. Obviously Brans - Dicke (BD) theory or its 
modifications have already found some attention as a driver of the present cosmic acceleration 
[331 El [351 EHl [371 HI] (see also [IS iH S!]). It had also been shown that BD theory can 
potentially generate sufficient acceleration in the matter dominated era even without any 
help from an exotic Q - field [32]. But this has problems with the required 'transition' from a 
decelerated to an accelerated phase. Amongst other nonminimally coupled theories, a dilatonic 
scalar field had also been considered as the driver of the present acceleration [45] . 

In most of the models the dark energy and dark matter components are considered to 
be non-interacting and are allowed to evolve independently. However, as the nature of these 
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components are not completely known, the interaction between them will indeed provide a 
more general framework to work in. Recently, Zimdahl and Pavon [46l |47j have shown that the 
interaction between dark energy and dark matter can be very useful in solving the coincidence 
problem ( see also ref [l8lH9| [50] ). Following this idea, we consider an interaction or 'transfer 
of energy' between the Brans - Dicke scalar field which is a geometrical field and the dark 
matter. The idea of using a 'transfer' of energy between matter and the nonminimally coupled 
field had been used earlier by Amendola [511 [52]. We do it specifically for a modified Brans - 
Dicke theory. The motivation for introducing this modification of Brans - Dicke theory is the 
following. In the presence of matter and a quintessence field, with or without an interaction 
between them, the evolution of net equation of state parameter w plays a crucial role in driving 
a late surge of accelerated expansion. But WMAP survey indicates that the time variation 
of w may be very severely restricted [53]- If the late acceleration is driven by an exchange of 
energy between matter and a geometrical field (j), the question of the variation of w would not 
arise. 

We write down the Brans - Dicke field equations in the so called Einstein frame. The field 
equations in this version look simpler and G becomes a constant. But one has to sacrifice 
the equivalence principle as the rest mass of a test particle becomes a function of the scalar 
field ^3j. So, the geodesic equation is no longer valid and the different physical quantities 
loose their significance. Nevertheless, the equations in this version of the theory enables us 
to identify the energy contributions from different components of matter. However, for final 
conclusions we go back to the original atomic units where we can talk about the features 
with confidence. We choose a particular form of the interaction and show that a constant 
BD parameter u can not give us the required fiip from a positive to a negative signature of 
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q in the matter dominated era. We attempt to sort out this problem using a modified form 
of BD theory where a; is a function of the scalar field (f) [55]. It has been pointed out by 
Bartolo and Pietroni [56] that a varying u theory can indeed explain the late time behaviour 
of the universe. By choosing a particular functional form of u;, we show that in the interacting 
scenario, one can obtain a scale factor 'a' in the original version ( i.e, in atomic units ) of 
the theory so that the deceleration parameter q has the desired property of a signature fiip 
without having to invoke any quintessence field in the model. We also calculate the statefinder 
pair {r,s}, recently introduced by Sahni et al [SZl [SH], for this model. The statefinder probes 
the expansion dynamics of the universe in terms of higher derivatives of the scale factor, i.e, 
d and a. These statefinder parameters along with the SNAP data can provide an excellent 
diagnostic for describing the properties of dark energy component in future. 

4.2.2 Field Equations and Solutions 

The field equations for a spatially flat Robertson - Walker spacetime in Brans - Dicke theory 



are 



3- = T + ii-3-T' (4-23) 

2^ + ^ = -^^^-^-2^^. (4.24) 

a a2 2 02 a0 ^ ' 

The field equations have been written with the assumption that at the present epoch the 
universe is filled with pressureless dust, i.e, Pm = 0. Here pm is the matter density of the 
universe, is the Brans - Dicke scalar field, a is the scale factor of the universe and uj is the 
BD parameter. An overhead dot represents a differentiation with respect to time t. 

The usual matter conservation equation has the form 
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p'm + SHpni = . 

But here we consider an interaction between dark matter and the geometrical scalar field and 
write down the matter conservation equation in the form 

Pm + 3Hp^ = Q , (4.25) 

such that the matter field grows or decays at the expense of the BD field. The matter itself is 
not conserved here and the nature of interaction is determined by the functional form of Q. 
We do not use the wave equation for the BD field here because if we treat equations ( ]4.23p . 
( 14.241) and (14.251) as independent equations, then the wave equation comes out automatically 
as a consequence of the Bianchi identity. It deserves mention that the wave equation will be 
modified to contain Q which will determine the rate of pumping energy from the BD field to 
matter or vice- versa. This interaction term Q is indeed a modification of Brans - Dicke theory. 
But this interaction does not demand any nonminimal coupling between matter and the scalar 
field and hence does not infringe the geodesic equation in anyway. In this interaction, the 
rest mass of a test particle is not modified but rather a "creation" of matter at the expense of 
the scalar field (or the reverse) takes place. In a sense, it has some similarity with the "C - 
field" of the steady state theory [59j. 

In the Brans - Dicke theory, the effective gravitational constant is given by G = which 
is indeed not a constant. Now, we effect a conformal transformation 

In the transformed version G becomes a constant. However, this transformation has some 
limitations which have been mentioned earlier. But the resulting field equations look more 



Acceleration of the Universe in Scalar - Tensor Theories 110 



tractable. Equations fl4.23p and (I4.24p in the new frame look like 



a a (2uj + 3) -n , _x 

2^ + - = , V ' , (4.27) 



and the matter conservation equation takes the form 

/>^ + 3-p^ = Q , (4.28) 
a 

where an overbar represents quantities in new frame and ip = Incj). The scale factor and the 
matter density in the present version are related to those in the original version as 

= (jyc? and pm = 4>^ pm ■ (4.29) 

Now, we choose the interaction Q of the form 

Q = -a Hpm , (4.30) 

where a is a positive constant. This negative Q indicates a transfer of energy from the dark- 
matter (DM) component to the geometrical field 0. 

Equation fl4.28p can be easily integrated with the help of equation fl4.30p to yield 

Pm = Po a^"""'^ , (4.31) 

where Pq is a constant of integration. 

Then, equations (14.261) and (14.271) alongwith equation (I4.3ip has a solution 

a = A?/(3+-) (4.32) 

where A is a constant given by 



Acceleration of the Universe in Scalar - Tensor Theories 111 



3-a V 2 

Some arbitrary constants of integration have been put equal to zero while arriving at equation 
(14.321) for the sake of simplicity. 

Using equations (I4.27P and (14.321) . one can easily arrive at the relation 

,2ct; + 3 -9 4a 1 , , 

If we consider a non- varying uj, equation fl4.33p will give rise to a simple power law evolution 
of (f). From equations (14. 29^ and (14.32^ . the scale factor a in atomic units will also have 
a power law evolution - an ever accelerating or an ever decelerating model contrary to our 
requirement. This is consistent with the exhaustive solutions in Brans - Dicke cosmology 
obtained by Gurevich et al. [50], where the dust solutions are all power law. This indicates 
that the chioce of constants of integration in equation (14.321) does not generically change the 
model. One way out of this problem is to consider a generalization of Brans - Dicke theory 
where the parameter is a function of the scalar field rather than a constant An 
evolving uj will be a contributory factor in determining the dynamics of the universe. 

We make a choice of u as, 

2CU + 3 a 



4 (3 + a)2(70-l)2 

Then, equation (14.331) can be integrated to yield 



(4.34) 



= (1 - 0ot)' , (4.35) 

00 being a positive constant. 

It deserves mention here that since g^o and ^oo are both equal to one, the time variable 
transforms as 
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df = (j)dt^ . 



This along with equation fl4.35p gives 



- 1 

00 



1 - ^20o(to - t) 



(4.36) 



which is a monotonically increasing function of t until t = Iq, beyond which the model really 
does not work. So one can use t itself as the new cosmic time in the original version of the 
theory without any loss of generality. So, for the sake of convenience, from now onwards, we 
write t in place of t. 

We transform the scale factor back to the original units by equation (14.291) . so that we are 
armed with the equivalence principle and can talk about the dynamics quite confidently. 
We have, 

a At2/(3+«) 

Also, the Hubble parameter and the deceleration parameter q in the original version comes 
out as, 

H = + , (4.38) 

3 + at 1- (pot ' ^ ' 

From equations (I4.37P and (I4.38P it is evident that at t — > both a and H blow up 
together giving a Big Rip. However, this rip has a different characteristic than that engineered 
by a normal phantom field. In the latter, goes to zero but pde goes to infinity at the rip. 
In the present case, however, there is no dark energy as such, and the scalar field is a part of 
geometry and hence it is difficult to recognize its contribution to the energy density. In the 
revised version, however, ^^^^ip^ is the contribution towards the stress tensor. It turns out 
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that at t — > this contribution remains quite finite. So the big rip is brought into being by 
the interaction, and not by a singularity in the stress tensor. 





1(a) 



1[b) 



Figure 4.1: Figure 4- i (a) and 4- i(b) shows the plot of q vs. t for different values of a . For figure 
4-l(ci) we choose a = 0.0000003 whereas for figure 4-^(b) we set the value as a = 0.0000001. 

The plot of q against t ( figure 4.1 ) reveals that the deceleration parameter indeed has a sign 
flip in the desired direction and indicates an early deceleration {q > 0) followed by a late time 
acceleration ( g < ) of the universe. Also, the nature of the curve is not crucially sensitive 
to the value of a chosen. 

Equation (14.271) clearly indicates that a/a is negative definite. So the signature fiip in 
Q = — ^2/^2 in the original Jordan frame must come from the time variation in via equation 



( I4.37p . Local astronomical experiments suggest that the present variation of G and hence that 
of has a very stringent upper bound. It is therefore imperative to check whether the present 
model is consistent with that bound. For figure 4.1, the value of the constant of integration 
00 is fixed at 0.3. With this value, and the age of the universe taken as ~ 15 Giga years, 
equation (14.351) yields 



10 per year. 
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which is consistent with the requirements of the local experiments [31] . The suffix indicates 
the present value. Also in this model, from equations (14.34^ and ( I4.37p . we get as u —>■ oo, 



from that in GR (a ~ t^/^). This is consistent with the notion that BD theory yields GR in 
the infinite u limit. 

4.2.3 Statefinder Parameters for the Model 

Recently Sahni et al. [571 [58] have introduced a pair of new cosmological parameters {r, s}, 
termed as "statefinder parameters". These parameters can effectively differentiate between 
different forms of dark energy and provide a simple diagnostic regarding whether a particular 
model fits into the basic observational data. These parameters are 



0^1 and a 



^2/(3+a)_ Therefore, for very small value of a (~ 10 ^), a is indistinguishable 



r = and s 



r-l 



3(9-1) • 



Accordingly, we find the statefinder parameters for the present model as 



r = 1 + 



3(j)lt^-3f3{l-(j)otf 2/3(1 - 0ot)3 + 20gt3 



(4.40) 



and 



s = 



3[0gt^ - P{1 - (j)otf][P + (1 - /3)0ot] + 2/3(1 - 0ot)3 + 20gt3 

3[/3 + (1 - /3)0ot] + (1 - mot} + Pil - 00^)2 - </)2t2]2 



(4.41) 




If we now plot r(s) for some small value of a ( a << 1 ), we find that the nature of the 
curve is similar to the one expected for scalar field quintessence models with equation of state 
parameter w in the range — 1 < w < [571 [58] . 
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Figure 4.2: Plot of r as a function of s for a = .0000001 . 
4.2.4 Discussion 

Thus we see that for a spatially flat universe {k = 0), we can construct a presently accelerating 
model in Brans - Dicke theory or more precisely in a generalized version of it ( = c<j(0) ) if 
one considers an interaction between dark-matter and the geometrical scalar field. It deserves 
mention that for this interaction, the Lagrangian should also be modified by the inclusion 
of an interference term between and L^- The salient feature of the model is that no dark 
energy component is required. The nature of the q vs. t curve is not crucially sensitive to small 
changes in the value of a, the parameter which determines the strength of the interaction; 
only the time of 'onset' of acceleration would change by small amounts with a. In revised 
unit, this interaction can be switched off by putting a = 0. However, in Jordan frame it is not 
possible to switch off the interaction with this particular choice. If we put a = 0, the scalar 
field itself becomes trivial. In this frame, the conservation equation ( equation fl4.28p along 
with equation (14.301) ) transforms to 

Pm . (4.42) 




p'm + 3Hp 



aH 



a 
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As one has both Hp^ and in the right hand side, the transfer of energy between matter 
and scalar field takes place both due to the expansion of the universe and the evolution of the 
scalar field. This equation shows that if 

(j) = constant a"^"/*^"^^-' 

the interaction vanishes. In this case, the transfer of energy due to Hpm and ^pm cancel each 
other. 

It is evident from equation (14.331) that if we consider the interaction between dark matter 
and the geometrical field of the form considered in equation (14.301) . a constant u will give rise 
to an ever accelerating or ever decelerating model and definitely we are not interested in that. 
So, the idea of varying uj is crucial here as it can very well serve the purpose of providing 
a signature flip in q in this interacting scenario. It deserves mention that the specific choice 
for the interaction in equation (I4.30p and the choice of a; = a;(0) in equation (I4.34p are 
taken so as to yield the desired result. This is indeed a toy model which simply shows that 
investigations regarding an interaction amongst matter and the nonminimally coupled scalar 
field is worthwhile. 

From equation (14.331) it is also evident that {2uj + 3) has to be positive definite, i.e, u has 
to pick up some positive value or at least uj should be greater than — | in order to sustain a 
consistent model. Also, the parameter uj does not have any stringent limit and thus it may be 
possible to adjust the value of uj to some higher value. Equation (14.341) indicates that if (p is 
very close to unity, which is consistent with the present value of G, uj can attain a high value 
at the present epoch. This is encouraging as it might be possible to obtain a model which 
exhibits early deceleration and late time acceleration even with a high value of u, compatible 
with the limit imposed on it by the solar system experiments |31j . 
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Also, from equation fl4.35l) it is evident that ^ < and | > 0. So, from equation (14.241) . 
which is of particular interest in studying the dynamics of the universe, we see that the last 
term 2^^ is negative and is the key factor in driving the present acceleration of the universe. 
This term basically provides the effective negative pressure and becomes dominant during 
the later stages of evolution and drives an accelerated expansion. We have also calculated 
the statefinder parameters for the model and show that the {r, s} pair mimics that of a 
quintessence model. For this however, the constant a should be given a very small value 
(~ 10^''). In this model we have considered a particular form of interaction which indeed is 
not unique, and some complicated kind of interaction may lead to more viable solutions of the 
various cosmological problems, particularly to a model which is not restricted in future. 

Although it is true that General Relativity (GR) is by far the best theory of gravity and 
the natural generalisation of BD theory to GR for large u limit is shown to be restricted in 
some sense |30j, still BD theory always seems to be ready to provide some useful clues to 
the solution of various cosmological problems. The solution to the "graceful exit" problem 
of inflation in terms of an 'extended inflation' scenario [211 |25] was first obtained in BD 
theory which provided hints towards the subsequent resolutions of the problem in GR. Once 
again here, BD theory in its own right could provide a model exhibiting the present cosmic 
acceleration without introducing any exotic dark energy component and since nothing definite 
is known about the source of this acceleration, this type of investigations may lead to some 
track along which viable solutions may finally be arrived at. However, this is a primitive 
model. This only shows that such investigations can be useful. It remains to be seen if the 
solution is an attractor, and whether the model is consistent with the structure formation. 
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Chapter 5 

Curvature Driven Acceleration 



1 oo 
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5.1 Curvature-driven Acceleration : a Utopia or a Re- 
ality? 

(Journal reference : S. Das, N. Banerjee and N. 
Dadhich, Class. Quantum Grav.,23, 4159 (2006).) 



[This work has been pubhshed in Classical and Quantum Gravity [Class. Quantum Grav., 
23, 4159 (2006)]. A minor error in the published version has been corrected in this version. 
The important conclusion that the model shows a signature flip in g at a finite past remains 
completely unaltered.] 
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5.1.1 Introduction 

The search for a dark energy component, the driver of the present accelerated expansion of 
the universe, has gathered a huge momentum because the alleged acceleration is now believed 
to be a certainty, courtesy the WMAP data [1]. As no single candidate enjoys a pronounced 
supremacy over the others as the dark energy component in terms of its being able to explain all 
the observational details as well as having a sound field theoretic support, any likely candidate 
deserves a careful scrutiny until a final unambiguous solution for the problem emerges. The 
cosmological constant A, a minimally coupled scalar field with a potential, Chaplygin gas or 
even a nonminimally coupled scalar field are amongst the most popular candidates ( see [2] for 
a comprehensive review). Recently an attempt in a slightly different direction is gaining more 
and more importance. This effort explores the possibility of whether geometry in its own right 
could serve the purpose of explaining the present accelerated expansion. The idea actually 
stems from the fact that higher order modifications of the Ricci curvature R, in the form of 

or R^iyR^'^ etc. in the Einstein - Hilbert action could generate an accelerated expansion 
in the very early universe [3]. As the curvature R is expected to fall off with the evolution, 
it is an obvious question if inverse powers of R in the action, which should become dominant 
during the later stages, could drive a late time acceleration. 

A substantial amount of work in this direction is already there in the literature. Capozziello 
et al. |1] introduced an action where R is replaced by i?" and showed that it leads to an 
accelerated expansion, i.e, a negative value for the deceleration parameter q for n = —1 and 
n = |. Carroll et al. [5J used a combination of R and and a conformally transformed version 
of theory, where the effect of the nonlinear contribution of the curvature is formally taken care 
of by a scalar field, could indeed generate a negative value for the deceleration parameter. 
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Vollick also used this 1/R term in the action [6] and the resulting field equations allowed 
an asymptotically exponential and hence accelerated expansion. The dynamical behaviour 
of gravity has been studied in detail by Carloni et.al |7]. A remarkable result obtained 
by Nojiri and Odinstov [8] shows that it may indeed be possible to attain an infiation at an 
early stage and also a late surge of accelerated expansion from the same set of field equations 
if the modified Lagrangian has the form L = R + i?™ + R~'^ where m and n are positive 
integers. However, the solutions obtained are piecewise, i.e, large and small values of the 
scalar curvature R, corresponding to early and late time behaviour of the model respectively, 
are treated separately. But this clearly hints towards a possibility that different modes of 
expansion at various stages of evolution could be accounted for by a curvature driven dynamics. 
Other interesting investigations such as that with an inverse sinh(i?) [9J or with InR terms 
[To] in the action are also there in the literature. 

The question of stability pjj and other problems notwithstanding, these investigations 
surely open up an interesting possibility for the search of dark energy in the non-linear contri- 
butions of the scalar curvature in the field equations. However, in most of these investigations 
so far mentioned, the present acceleration comes either as an asymptotic solution of the field 
equations in the large cosmic time limit, or even as a permanent feature of the dynamics of 
the universe. But both the theoretical demand [12] as well as observations \i3\ ( see also py 
) clearly indicate that the universe entered into its accelerated phase of expansion only very 
recently and had been decelerating for the major part of its evolution. So the deceleration 
parameter q must have a signature fiip from a positive to a negative value only in a recent 
past. 

In the present work, we write down the field equations for a general Lagrangian f{R) and 
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,4 



investigate the behaviour of the model for two specific choices of f{R), namely f{R) = R — ^ 
and f{R) = . 

Although the field equations, a set of fourth order differential equations for the scale 
factor a, could not be completely solved analytically, the evolution of the 'acceleration' of the 
universe could indeed be studied at one go, i.e, without having to resort to a piecewise solution. 
The results obtained are encouraging, both the examples show smooth transitions from the 
decelerated to the accelerated phase. In this work we virtually assume nothing regarding the 
relative strengths of different terms and let them compete in their own way, and still obtain the 
desired transition in the signature of the deceleration parameter q. This definitely provides a 
very strong support for the host of investigations on curvature driven acceleration, particularly 
those quoted in [5, 6, 7 and 8]. 

In the evolution equation, q is expressed as a function of if, the Hubble parameter. This 
enables one to write an equation with only q to solve for; as the only other variable remains 
is H which becomes the argument. This method appears to be extremely useful, although it 
finds hardly any application in the literature. The only example noted by us is the one by 
Carroll et al. [2], which, however, describes the nature only in an asymptotic limit. 

In the next section the model with two examples are described and in the last section we 
include some discussion. 

5.1.2 Curvature Driven Acceleration 

The relevant action is 




(5.1) 
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where the usual Einstein - Hilbert action is generahzed by replacing R with f{R), which is an 
analytic function of R, and Lm is the Lagrangian for all the matter fields. A variation of this 
action with respect to the metric yields the field equations as 

Cr R^i/ '^Rq^v + , (5-2) 

where the choice of units SttG = 1 has been made. T^i,^ represents the contribution from 
matter fields scaled by a factor of jjj^ and T^^,^ denotes that from the curvature to the effective 
stress energy tensor. T^^'^ is actually given as 



\9tiu{f{R) - Rf'{R)) + - atiuQap) 



(5.3) 



.2^ 

The prime indicates differentiation with respect to Ricci scalar R. It deserves mention that we 
use a variation of (15. ip w.r.t. the metric tensor as in Einstein - Hilbert variational principle and 
not a Palatini variation where A is varied w.r.t. both the metric and the affine connections. 
As the actual focus of the work is to scrutinize the role of geometry alone in driving an 
acceleration in the later stages, we shall work without any matter content, i.e, Lm = leading 
to T^jy^^ = 0. So for a spatially flat Robertson - Walker spacetime, where 

ds^ = dt^ - a^{t)[dr'^ + r'^de'^ + sin^ ^rf^^], (5.4) 

the field equations (15. 2p take the form ( see [3] ) 



1 

a- ~ r 



hf-Rf')-3^Rf" 
2 a 



(5.5) 



2- + ^ = -^[2-i?r + Rf" + R'f" -kf- Rf')] . (5.6) 

Here a is the scale factor and an overhead dot indicates differentiation w.r.t. the cosmic time 
t. If f{R) = R, the equation (15.20 and hence (15.50 and (15. 6p take the usual form of vacuum 
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Einstein field equations. It should be noted that the Ricci scalar R is given by 



R = -6 



• 2 ■ 

a a 

- + - 

a 



(5.7) 



and already involves a second order time derivative of a. As equation fl5.6l) contains i?, one 
actually has a system of fourth order differential equations. 

It deserves mention at this stage that if i? is a constant, then whatever form of f{R) 
is chosen except f\R) = R, equations (15.51) and (15.61) represent a vacuum universe with a 
cosmological constant and hence yield a deSitter solution, i.e, an ever accelerating universe. 
Evidently we are not interseted in that, we are rather in search of a model which clearly shows 
a transition from a decelerated to an accelerated phase of expansion of the universe. As we 
are looking for a curvature driven acceleration at late time, and the curvature is expected to 
fall off with the evolution, we shall take a form of f{R) which has a sector growing with the 
fall of R. We work out two examples where indeed the primary purpose is served. 



(i) f{R) = R-^. 



In the first example, we take 



fiR) = R 



R 



where /i is a constant. Indeed fi has a dimension, that of R^, i.e, that of (time)^^. This is 
exactly the form used by Carroll et al. [5] and Vollick [6]. Using the expression (15.81) in a 
combination of the field equations (15.51) and (15.61) . one can easily arrive at the equation 

1 



2H 



(i?2 + 



(5.9) 
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where H = ^, is the Hubble parameter. As both R and H are functions of a and its derivatives, 
equation (I5.9P looks set for yielding the solution for the scale factor. But it involves fourth 
order derivatives of a ( i? already contains d ) and is highly nonlinear. This makes it difficult 
to obtain a completely analytic solution for a. As opposed to the earlier investigations where 
either a piecewise or an asymptotic solution was studied, we adopt the following strategy. The 
point of interest is the evolution of the deceleration parameter 

aa H 

q = -^ = ——-l. (5.10) 



a 



So we translate equation (15.91) into the evolution equation for q using equation (I5.10p and obtain 

l^'j^) - 3/i^(^ - 6l^'H\q + if + 3f,^H\q + 1) 

+ 3^^H^j^^^q + 36H%q - ifiq + 1) = 0. (5.11) 

This equation, although still highly nonlinear, is a second order equation in q. But the problem 
is that both q and H are functions of time and cannot be solved for with the help of a single 
equation. However, they are not independent and are connected by equation ( ]5.10p . So we 
replace time derivatives by derivatives w.r.t. H using equation (I5.10p and write (15. lip as 

^(g2-l)ifV^-^if2(g+2)/-lif(g-l)(4g+7)gt-(2g+l)(g - lf + H\q-l)' = 0. (5.12) 

Here for the sake of simplicity /i^ is chosen to be 12 ( in proper units ), and a dagger represents 
a differentiation w.r.t. the Hubble parameter H. As is a measure of the age of the universe 
and if is a monotonically decreasing function of the cosmic time, equation (I5.12p can now be 
used as the evolution equation for q. The equation appears to be hopelessly nonlinear to give 
an analytic solution but if one provides two initial conditions, for q and g^, for some value of 
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H, a numerical solution is definitely on cards. We choose units so that i^o, the present value 
of H, is unity and pick up sets of values for q and for H = 1 { i.e, the present values ) from 
observationally consistent region [15] and plot q versus H numerically in figure 5.1. As the 
inverse of H is the estimate for the cosmic age, 'future' is given hj H < 1 and past by if > 1. 
The plots speak for themselves. One has the desired feature of a negative g at if = 1 and 
it comes to this negative phase only in the recent past. An important point to note here is 
that neither the nature of the plots, nor the values of H at which the transition takes place, 
crucially depends on the choice of initial conditions, so the model is reasonably stable. 




0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 

>H >H 



Figure 5.1: Figure 5.1(a) and 5.1(h) shows the plot of q vs. H for f{R) = R— ^ for different 
initial conditions. For figure 5.1(a) we choose the initial conditions as q[l] = —0.5, q'[l] = 1.2 
whereas for figure 5.1(b) we set the initial conditions as q[l] = —0.5, q'[l] = 1.0. 

(ii) f{R) = e-f . 



In this choice, the function / is monotonically increasing with t as i? is decreasing with t. 
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The field equations (15. 5p and (15.61) have the form 



• 2 

a a 
2- + — = -6 
a 



2* 6' 



2^ 6' 



1 d 
12a' 



-R 



Id- 1 •■ 1 • 9' 

R R+ R^ 

18 a 36 216 



(5.13) 
(5.14) 



From these two equations it is easy to write 



2H = -R- —R^ - -HR. 
6 36 6 



(5.15) 



Following the same method as before, the evolution of g as a function of H can be written as 



H\q + l)gtt + [H^ _ (g + l)H^]qt'^ + [2{q + 1)H^ + 6qH^ + 3H^ - 4H^{q'^ - l)]gt 

+ 6H\'^ + 2H'^q - 4H\q^ - l){q - 1) - 8H^ + 2 = . (5.16) 

With similar initial conditions for q and at if = 1, the plot of q versus H ( figure 5.2 ) 
shows features similar to the previous example, the deceleration parameter q has a signature 
change, from a positive to a negative phase in the recent past { H > 1 ). This example has 
an additional feature that the universe re-enters a decelerated phase of expansion again in a 
near future [H < 1). In this case also, a small change in initial conditions hardly has any 
perceptible change in the graphs. 

In this curvature singularity in a finite future is indicated. As g = —-^j^ and 

H remains finite, ^ and hence the curvature (via equation (15. 7p ) has a singularity in a finite 
future. This is consistent with [5] which indicates that a curvature quintessence may end 
up with three possibilities - an asymptotic de Sitter, a power law infiation or a curvature 
singularity in a finite future. The present case corresponds to the third possibility. The 
curves however show quite clearly that this singularity is not a 'Big Rip' type, where due to 
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continuous vigorous acceleration both H and a blow up in some finite future. Here the model 
clearly enters into a decelerating phase close to H = 0.8, as shown by the figure. 

From equation (15.161) . one can also conclude that for very high value of H (i.e, when the 
age of the universe was very small), q —>■ —1, which gives an early inflation. 




o.e o.s 1 i.£ 1.4 i.e o.f 5,5 j, j,,^ j,,^ 

" H 



Fig 2(a) Fig 2(b) 

Figure 5.2: Figure 5.2(a) and 5.2(h) shows the plot of q vs. H for f{R) = for different 
initial conditions. Here also for figure 5.2(a) and 5.2(b) we set the initial conditions as q[l] = 
—0.5, q'[l] = 1.0 and q[l] = —0.5, q'[l] = 1.2 respectively. 

As the plots provide a sufficient data set, attempts could be made to find the closest 
analytical expression for q = q{H). These expressions are found to be polynomials. For 
example, a very close analytical expression for figure 5.2(b), within the accuracy of plots, is 
given as 

q = 47.95H^ - 335.73i/^ + 991H^ - 1586.90i/^ + 1459.20i/^ - 729.73i/ + 153.74. (5.17) 

This expression holds only when H is reasonably close to one, and has nothing to do with 
other ranges of values of H. 
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5.1.3 Discussion 

The present work indicates that by asking the question whether geometry in its own right 
can lead to the late surge of accelerated expansion, some feats can surely be achieved. Both 
the examples considered here indicate that one can build up models which start accelerating 
at the later stage of evolution and thus allow all the past glories of the decelerated model 
like nucleosynthesis or structure formation to remain intact. An added bonus of the second 
example is that the universe re-enters a decelerated phase in near future and the 'phantom 
menace' is avoided - the universe does not have to have a singularity of infinite volume and 
infinite rate of expansion in a 'finite' future. 

It is of course true that a lot of other criteria have to be satisfied before one makes a final 
choice, and we are nowhere near that. Already there is a criticism of gravity that it is 
unsuitable for local astrophysics because of problems regarding stability However, it was 
pointed out by Nojiri and Odinstov [S] that a polynomial may save the situation ( see also 
reference [16] ). Our second example is exponential in R, i.e, a series of positive powers in R, 
and hence could well satisfy the criterion of stability. As already pointed out, although the 

4 

choice of f{R) = R—^ is already there in the literature and served the purpose in a restricted 
sense than it does in the present work, the choice of f{R) = has hardly any mention in 
the literature. The Lagrangian f{R) = contains a cosmological constant as f{R) ~ 1 — 
for small R. So indeed one expects that it gives an accelerated expansion. But the interesting 
feature is that the same model gives an early inflation followed by a decelerated expansion, 
then an accelerated expansion around the present epoch and a decelerated phase once again 
in near future. 

It should also be noted that the present toy model deals with a vacuum universe and one 
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has to either put in matter, or derive the relevant matter at the right epoch from the curvature 
itself. Some efforts towards this have already begun [17]. On the whole, there are reasons 
to be optimistic about a curvature-driven acceleration which might become more and more 
important in view of the fact that WMAP data could indicate a very strong constraint on the 
variation of the equation of state parameter w [18] . 
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